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Abreact 


This  report  describe  a compute  program  that  simulates  the  growth  and  dtclina  of 
bacterial  population*  In  food*  during  Pta  vartou*  opantion*  of  food-praparadon  and 
praoorinp  By  means  of  this  program  an  existing  or  propoaad  food  pnparatlon  procaw 
can  ba  ttudlad,  ha  wfaty  awawad  and  tha  critical  oparatiom  identified.  Also  hypothetical 
systems  for  monitoring  tha  bacterial  growth  can  ba  triad  h the  simulation.  Tha  basic 
hypothad*  of  tha  aimulatlon  is  that  bacterial  growth  dapands  on  tima  and  temparature 
in  a food.  Detailed  Instructions  for  die  use  of  tha  program  are  given,  and  an  example 
is  set  up  and  solved. 


I.  Introduction 


In  recent  years  the  U.  S.  Army  has  undertaken  a major  effort  to  improve  its 
food-service  system,  especially  the  garrison  food-service  system  which  affects  the  morale 
of  enlisted  men  so  directly.  The  existing  system  has  been  studied  extensively  and  various 
proposals  made  for  its  moderization.  Some  of  the  suggested  changes  in  the  system  are 
new  enough  so  there  is  no  backlog  of  civilian  experience  as  a guide  to  their  bacteriological 
safety.  This  situation  was  perceived  by  members  of  the  Operations  Research  and  Systems 
Analysis  Group  and  the  Food  Laboratory  at  the  U.  S.  Army  Natick  Laboratories,  and 
a method  was  sought  that  would  permit  rational  assessment  of  the  bacteriological  safety 
of  these  processes.  The  present  work  is  an  attempt  by  the  author  to  provide  such  a 
method. 

This  report  describes  a computer  program  that  simulates  the  growth  and  decline  of 
bacterial  populations  in  foods  during  the  various  operations  of  food-preparation,  i.e., 
storing,  mixing,  cooking,  cooling,  etc.  It  is  useful  for  at  least  two  purposes.  First,  it 
can  simulate  in  existing  or  proposed  food-preparation  process  and  tell  whether  the  process 
is  safe  or  not.  Doing  this  permits  the  user  to  see  where-abouts  in  the  process  the  most 
dangerous  situation  occurs  and,  sometimes,  what  can  be  done  to  remedy  it.  Second  it 
provides  a means  for  testing  the  effectiveness  of  systems  for  monitoring  the  bacteriological 
safety  of  food-preparation  processes. 

We  must  emphasize  that  in  order  to  simulate  a food-preparation  process,  that  process 
must  be  very  well-defined,  much  more  so,  for  example,  than  is  done  by  the  recipe  for 
the  food.  This  means  that  the  user  may  have  to  gather  a lot  of  information  before  running 
a simulation.  This  is  of  course  inconvenient,  but  it  may  have  benefits,  too.  A user 
of  the  simulation  who  is  investigating  a potentially  dangerous  situation  in  preparing  a 
certain  food  may  find  that  merely  defining  the  process  carefully  is  enough  to  let  him 
see  where  the  trouble  may  occur  and  how  serious  it  is. 

The  computer  program  is  written  in  the  FORTRAN  V language  rather  than  in  one 
of  the  languages  specifically  invented  for  simulation  purposes,  such  as  GPSS,  SIMSCRIPT, 
SIMULA  or  GASP.  This  choice  was  dictated  partly  by  the  unavailability  of  some  of 
these  compilers  for  the  UNIVAC  1106,  (the  Natick  Laboratories  computing  system)  as 
well  as  the  author's  background,  the  pressure  of  time  and  the  nature  of  the  simulation. 


The  principal  hypothesis  on  which  the  model  is  based  is  that  growth  (or  death)  of 
bacteria  in  a food  during  storage  depends  on  time  and  temperature.  This  dependence 
is  modelled  in  the  computer  program  by  general  formulas  containing  six  constants  for 
each  stored  food.  The  values  assigned  to  these  constants  specify  the  growth  and  death 
behavior  of  the  bacteria  in  that  food.  The  model  can  easily  handle  the  situation  where 
the  bacterial  growth  is  different  in  the  various  foods,  provided  the  differences  are  ones 
that  can  be  accurately  modelled  by  changes  in  the  constant  values. 

At  first  glance  it  appears  that  die  principal  hypothesis  is  unduly  restrictive  in  two 
respects.  First,  only  one  organism  can  be  studied.  Second,  bacterial  growth  depends 
on  many  factors  other  than  time  and  temperature,  such  as  moisture  content,  ph,  salinity, 
etc.  In  fact  neither  of  these  apparent  difficulties  is  too  serious,  for  they  can  usually 
be  avoided  by  repeated  running  of  the  program  using  different  values  of  the  constants 
for  the  various  foods.  For  example  the  growth  of  two  different  organisms  in  a food 
process  can  be  studied  by  making  two  computer  runs,  assigning  different  values  to  the 
sets  of  constants  that  prescribe  the  growth  properties  of  the  two  organisms.  This  procedure 
will  usually  be  adequate  but  may  fail  if  the  two  organisms  interact  significantly  with  each 
other.  Again,  the  effect  of  two  different  moisture  levels  can  be  studied  provided  the 
growth  of  the  organism  at  different  moisture  levels  can  be  accurately  described  by  changes 
in  the  values  of  the  six  constants  that  define  growth  behavior.  This  will  often  give 
accurate-enough  results  but  could  be  misleading  if  growth  is  sensitive  to  moisture  and 
if  the  moisture  in  one  or  more  foods  changes  radically  with  time. 

The  general  structure  of  the  model  is  described  in  Section  2 and  a brief  outline 
of  the  computer  program  is  given  in  Section  3.  Section  4 describes  an  example  of  the 
use  of  the  program  and  the  preparation  of  the  data.  The  results  of  this  example  and 
some  further  aspects  of  the  model  and  program  are  discussed  in  Section  5.  Appendix  I 
contains  the  derivation  of  the  principal  formulas.  In  Appendix  II  rhere  is  a detailed 
description  of  the  input  data,  keyed  to  the  example  described  in  Section  4.  The  printout 
of  the  results  of  the  example  forms  Appendix  III,  and  Appendix  IV  contains  a listing 
of  the  FORTRAN  V program.  The  set-up  of  executive  control  cards  to  be  used  when 
running  the  program  on  the  Natick  Laboratories  UNI  VAC  1106  System  is  described  in 
Appendix  V. 
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2.  General  Structure  of  the  Model 


The  main  general  types  of  entities  that  occur  in  the  mode!  are  batches  of  food, 
having  one  main,  and  several  auxiliary,  attributes  of  interest,  and  activities  performed  on 
these  batches  of  food  that  cause  changes  in  the  attributes.  The  independent  variable, 
time,  also  plays  an  important  part.  We  shall  describe  the  roles  of  these  entities  in  increasing 
detail  in  the  following  paragraphs. 

The  independent  variable,  time,  is  defined  by  ihe  discrete  (integer-valued)  quantity  I. 
The  simulation  starts  at  I = 1 and  proceeds  until  I attains  a value  not  exceeding  IMAX. 

Each  batch  (packet)  of  food  is  assigned  an  index,  J,  and  has,  as  its  main  attribute, 
bacteria  count  and,  as  auxiliary  attributes,  temperature  and  weight.  The  bacteria  count 
is  specified  by  the  quantity  V(J,  I),  which  is  the  natural  logarithm  of  the  density  of 
organisms  (e.g.  number  of  organisms  per  gram)  in  the  J-th  food  at  the  l-th  time-step. 
The  temperature  of  the  J-th  batch  of  food  at  time-step  I is  Y(J,  I),  and  the  weight  of 
the  batch  is  W(J),  which  is  assumed  not  to  depend  on  the  time. 

The  two  principal  types  of  activities  that  are  performed  are  storage  and  food -hand ling. 
Each  storage  area  is  assigned  an  index,  KS,  and  a temperature,  YST(KS),  which  is  the 
ambient  temperature  of  all  food  batches  in  that  storage  area.  Every  food  is  in  some 
storage  area  while  it  is  in  the  system.  It  is  important  to  notice  that  we  regard  cooking 
merely  as  storage  at  high  temperature,  Food-handling  activities  consist  of  two  operations, 
namely  splitting  and  blending  batches  of  food.  In  the  splitting  operation,  a single  batch 
of  food  is  subdivided  into  several  new  foods,  and  in  the  blending  operation  several  batches 
of  food  are  mixed  together  to  give  one  new  food-batch.  Splitting  and  blending  are  each 
other's  mirror  images  with  respect  to  time  and  are  assumed  to  occur  instantaneously,  i.e. 
in  much  less  than  one  time-step. 

Foods  may  be  in  the  syster  at  the  start  of  simulation  or  brought  in  at  a later  time, 
and  in  either  case  the  initial  values  of  the  attributes  must  be  prescribed  by  the  user. 
Foods  may  also  be  created  or  born  by  means  of  splits  or  blends  of  foods  previously 
in  the  system,  in  which  case  the  values  of  the  attributes  are  derived  from  those  of  the 
progeritor  foods  by  the  formulas  described  in  Appendix  I.B.  Foods  may  leave  the  system 
by  giving  birth  to  other  foods  at  a split  or  blend,  or  they  may  be  served,  or  they  may 
remain  in  the  system  (i.o.  in  some  storage)  until  the  end  of  the  simulation. 


3 


Bacteria  counts  may  change  due  to  either  natural  growth  (i.e.  multiplication)  and 
death  or  contamination.  Natural  growth  or  death  is  assumed  to  occur  as  time  passes 
during  storage.  The  ratBS  of  change  per  time-step,  A(l),  depends  on  the  present  bacteria 
count  and  temperature  of  the  food.  The  temperature  of  the  food  is  first  calculated  from 
its  initial  temperature,  ambient  temperature  and  effective  diffusion  coefficient.  This 
temperature  is  then  employed  in  general  formulas  for  A(l),  and  a simple,  approximate, 
numerical  integration  of  A(l)  is  used  to  find  the  bacteria  count,  V(J,  I).  The  quantity 
A(l)  describes  both  the  growth  rate  (when  A(l)  > 0)  and  the  death  rate  (when  A(l)  < O) 
of  the  bacteria.  The  formulas  for  calculating  A(l)  contain  a number  of  constants  that 
the  user  must  provide  in  order  to  model  the  growth  behavior  of  the  organism  being  studied. 
The  details  of  the  formulas  for  calculating  natural  growth  and  death  of  organisms  are 
given  in  Appendix  I. 

The  other  source  of  change  in  bacteria  counts  is  the  contamination  introduced  into 
the  food.  This  may  occur  gradually  in  a stored  food  or  instantaneously  at  a split  or 
blend  and  is  modelled  by  generating  a positive  random  number  and  adding  it  to  the  value 
of  V.  For  a food  created  at  a split  or  blend  this  is  done  after  the  values  of  V implied 
by  the  progenitor  foods  have  been  calculated. 

In  addition  to  these  activities  the  model  also  simulates  two  systems  for  monitoring 
the  number  of  organisms.  The  first  is  a record-keeping  system  that  keeps  track  of  the 
lengths  of  time  each  food  has  been  exposed  to  its  ambient  temperatures,  estimates  the 
bacteria  count  at  each  handling,  compares  the  estimate  with  a pre-assigned  threshold  and 
gives  a signal  if  the  threshold  value  is  exceeded.  The  second  system  is  intended  to  model 
the  domestic,  or  housewife,  method.  At  each  handling  the  threshold  bacteria  count  is 
compared  with  the  counts  of  all  the  foods  directly  involved  in  the  handling  or  located 
in  the  same  storage  as  any  directly  - involved  food.  Again  a signal  is  given  if  the  threshold 
value  is  exceeded. 

The  system  makes  allowance  for  errors  of  various  types  of  means  of  random  numbers, 
usually  obeying  a Gaussian  (Normal)  distribution.  The  principal  kinds  of  errors  in  the 
model  are  scheduling  errors,  errors  in  the  attribute  values  of  foods  initially  in  the  system 
or  entering  it  from  outside,  errors  in  storage  temperatures,  random  losses  of  weight  (spillage) 
and  increases  in  contamination  during  handling  (previously  mentioned),  and  variation  in 
sensitivity  of  the  cook  in  the  domestic  monitor  system. 
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Finally  the  model  also  permits  changes  in  the  ambient  temperature  of  each  food. 
Up  to  five  changes  to  new  ambient  temperatures  at  arbitrary  times  may  be  made  in  each 
food,  and  random  errors  may  occur  in  the  times  of  these  changes. 

3.  Brief  Sketch  of  Computer  Program 

The  computer  program  which  carries  out  the  simulation  consists  of  three  major  parts 
and  two  important  subroutines.  The  first  major  part  is  the  reading  and  recording  of 
the  input  data.  For  each  food  the  input  data  consists  of  such  information  as  scheduling, 
storage  area,  initial  attribute  values,  constants  defining  growth  rates  and  data  about  random 
errors.  Also  for  each  split  or  blend  information  is  given  about  which  foods  are  involved 
in  the  operation,  random  errors  in  Its  time  of  occurence  and  the  handling  contamination. 
In  addition  certain  general  data  such  as  the  number  of  repetitions  of  the  simulation  and 
the  temperatures  of  the  storage  areas  are  also  read. 

The  second  main  portion  consists  of  deriving  necessary  scheduling  information  from 
the  input  data  and  introducing  the  random  errors  into  that  scheduling  and  into  the  values 
of  V,  Y and  W fcr  the  foods  initially  in  the  system  or  entering  it  from  outside.  For 
each  repetition  of  the  simulation  this  results  in  a final  schedule  and  set  of  input  values 
for  the  attributes  possibly  somewhat  different  from  the  schedule  and  values  originally 
intended.  The  linal  schedule  and  values  are  printed  at  the  start  of  each  repetition  of 
the  simulation. 

The  third  main  part  of  the  program  is  the  actual  calculation  of  V,  Y and  W during 
both  storage  and  handling.  This  is  largely  done  by  means  of  the  two  important  subroutines, 
BACGRO  and  SPLEND.  BACGRO  does  the  entire  calculation  and  printing  out  of  Y, 
V,  A and  the  ambient  temperature  at  each  time-step  during  the  storage  of  each  food. 
SPLEND  carries  out  the  calculations  of  V,  Y and  W for  the  new  foods  created  at  a split 
or  blend  and  incorporates  the  random-handling  contamination  and  spillage  losses  into  the 
values  of  V and  W,  respectively. 

The  instructions  that  simulate  the  two  monitoring  activities  are  spread  throughout 
the  program  and  can  be  identified  by  the  annotations  MON R EC  for  the  record-keeping 
system  and  MON  DOM  for  the  domestic  system.  The  short  subroutine  SNIFF  is  part 
of  the  domestic  system  and  carries  out  the  determination  of  whether  the  cooks  sense 
that  a food  is  "bad"  or  not. 


5 


4.  Ex  am  pi* 


We  shall  describe  a simulation  of  the  preparation  of  chicken  pot  pie,  following  a 
recipe  similar  to  that  in  the  Army  CFPF  production  guide  for  pre-cooked,  frozen  and 
then  reheated  chicken  pot  pie.  This  is  a moderately  complicated  process  and  illustrates 
a numbe-  of  the  features  of  the  program  as  well  as  some  of  the  problems  that  arise  in 
translating  a recipe  into  data  that  the  program  can  use. 

The  first  sttp  in  modelling  this  process  is  to  study  the  recipe,  list  the  ingredients 
and  try  to  get  a clear  idea  of  what  operations  are  performed  and  when  they  occur  in 
the  process.  Frequently  the  recipe  is  not  clear  about  some  of  the  details  of  the  process, 
and  one  must  either  seek  more  information  or  try  to  make  sensible  guesses.  In  the  present 
example  we  start  with  a fairly  detailed  list  of  ingredients: 

60  lbs  of  chicken 

2 lbs  of  sauce  (roux) 

23  lbs  of  stock  and  seasonings 

8 lbs  of  chilled  vegetables 

6 lbs  of  frozen  peas 

2 lbs  of  slurry,  made  of  starch  and  water 
5 lbs  of  topping  dough 

Now  it  is  clear  that  in  listing  the  ingredients  in  this  way,  we  have  already  made  some 
decisions  about  the  amount  of  detail  that  we  want  our  model  to  have.  For  example 
we  have  decided  not  to  bother  with  the  details  of  how  the  slurry  is  made  from  starch 
and  water  and  how  much  seasoning  of  various  kinds  goes  into  the  stock.  Typically  in 
starting  a model  a number  of  such  decisions  have  to  be  made,  in  which  the  gains  from 
added  information  must  be  weighed  against  the  cost  of  greater  complexity. 
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Symbol  for  a split,  occurring  at  1*5,  in  which  food  1 is  split  into  foods  3,  7 and 
8.  The  split  number  is  NS=2. 
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Symbol  for  a blend,  occurring  at  1=12,  in  which  foods  4,  5 and  7 are  combined 
to  give  food  8.  The  blend  number  is  NB=3.  Notice  that  food  6 does  not  take  part 
in  the  blend. 


Symbol  for  serving  of  food  number  9 at  time  1=83. 


Figure  1:  Diagrammatic  Symbols  for  Splits,  Blends  and  Servings 
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The  activities  of  the  recipe  can  be  described  as  follows:  in  the  central  preparation 
area  cook  the  chicken,  remove  the  chicken  from  the  juices  (stock)  and  let  it  cool.  Combine 
this  stock  with  the  sauce,  chilled  vegetables  and  stock  and  seasoning  prepared  in  advance. 
Heat  the  mixture,  add  frozen  peas  then  heat  again  and  add  slurry  and  the  cooked  chicken. 
Heat  and  add  topping,  neat  again,  then  allov/  it  to  cool.  The  pios  are  frozen  and  transported 
to  the  serving  area,  where  they  are  heated  and  served. 

This  description  is  incomplete  since  it  includes  no  quantitative  information  about 
storage  temperatures  or  times.  After  further  investigation,  a description  was  reached  which 
is  best  expressed  by  the  diagram  shown  in  Figure  2.  This  diagram  uses  the  symbols  for 
splits,  blends  and  servings  described  in  Figure  1.  In  this  diagram  the  horizontal  scale 
is  time-steps  of  five  minutes  each. 

Thirteen  foods  in  all  occur  during  the  preparation,  of  which  five  are  present  at  the 
initial  time,  namely  foods  numbered  1 (chicken),  2 (sauce),  3 (stock  and  seasoning),  5 
(chilled  vegetables),  7 (frozen  peas).  Two  foods,  numbered  9 (slurry)  and  12  (topping) 
are  brought  in  at  time  steps  I = 25  and  I = 35,  respectively,  i.e.  about  125  and  175 
minutes  after  the  start  of  the  process.  Six  foods  are  created  during  the  process  and  are 
numbered  4 (extra  stock),  10  (de-boned  chicken),  13  (final  chicken  pot  pie),  6 (alpha), 
8 (beta)  and  11  (gamma).  The  last  three  are  simply  mixtures  created  at  intermediate 
stages  during  the  food  process. 

Each  food  is  represented  by  a horizontal  line  in  Figure  2,  extending  over  the  time-steps 
at  which  that  food  exists.  Each  line  ends  at  a split  or  blend  or  serving  symbol.  The 
storage  unit  number  for  that  food  is  written  above  its  line,  and  the  temperatures  of  the 
various  storage  units  are  shown  in  the  upper  right  comer  of  Figure  2.  Changes  in  the 
ambient  temperature  are  shown  by  writing  the  new  ambient  temperature  to  the  right  of 
a short,  vertical,  dashed  line.  For  example,  we  see  that  food  13,  the  final  chicken  pot 
pie  is  first  cooled  at  room  temperature,  since  it  is  in  storage  unit  3,  from  time  I = 39 
to  I = 45.  It  is  put  in  the  freezer  (which  has  ambient  temperature  —18°)  from  I = 45 
to  I = 91,  then  reheated  at  temperature  175°  from  I = 91  to  97  and  served. 


8 


4 


ST.  UNIT  TEMPERATURES 
ST  1 is  at  -18°  ' 

ST  2 is  at  5°C 
ST  3 is  at  20°C 
ST  4 is  at  100°C 
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Beta 

9 

Slurry 

10 

De-Boned  Chicken 
11 

Gamma 

12 

Topping 

13 

Chicken  Pot  Pie 


In  the  author's  experience  it  is  usually  very  helpful,  almost  necessary,  to  prepare 
a diagram  like  Figure  2 if  the  process  being  studied  is  at  all  complicated.  With  such 
a diagram  it  is  fairly  straightforward  to  prepare  the  input  data  although  some  additional 
information,  such  as  the  initial  attribute  values,  has  to  be  added  to  that  of  the  diagram. 
In  the  present  example  the  initial  temperature,  V,  bacteria  count,  V,  and  weight,  W,  are 
listed  in  Table  I for  each  food  that  is  not  created  from  existing  foods  during  the  process. 

Appendix  II  contains  a complete  description  of  the  input  data  for  this  simulation. 
Five  repetitions  of  the  entire  simulation  are  run,  with  different  random  numbers  at  each 
repetition.  The  output  from  one  of  the  repeti*:ons  is  listed  in  Appendix  III.  Graphs 
of  Y and  V as  functions  of  time  for  certain  of  the  foods  are  shown  in  Figure  3. 

5.  Discussion 

The  principal  results  of  the  chicken  pot  pie  simulation  are  listed  in  Appendix  III, 
arising  from  the  input  data  of  Appendix  II.  Some  of  the  numbers  used  in  this  example, 
(particularly  those  connected  with  the  action  of  the  monitor  systems)  are  not  realistic 
because  they  were  chosen  to  demonstrate  certain  aspects  of  the  program.  For  example, 
the  constants  P(J,  K),  describing  the  growth  of  the  organisms,  are  very  roughly  those 
of  E.  coli  at  five-minute  time-steps.  It  is  possible  that  we  would  not  regard  food  as 
being  truly  spoiled  if  it  contained  1200  organisms  per  gram,  which  is  what  we  mean  by 
choosing  TRT  = 7.150.  The  choice  TRT  * 7.150  was  made  to  ensure  that  we  obtained 
some  cases  that  we  could  regard  as  truly  spoiled. 

Examination  of  the  results  of  the  simulation  show  that  the  final  food,  number  13, 
chicken  pot  pie,  has  very  low  bacteria  counts,  much  less  than  one  per  gram.  The  highest 
bacteria  counts  were  always  obtained  in  the  chilled  vegetables,  food  number  5,  mostly 
because  these  were  brought  into  the  system  with  a high  count  (V  = 7.0  initially).  In 
a sense,  therefore,  this  example  is  a rather  uninteresting  one  because  the  food  that  is 
served  is  very  safe. 

A possible  uje  of  this  simulation  is  to  investigate  how  "stable"  this  process  is  by 
seeing  how  much  the  results  change  when  we  alter  some  of  the  input  data.  For  example, 
we  might  greatly  reduce  the  final  reheating  by  setting  the  last  storage  temperature  for 
food  number  13  at  50°C  instead  of  175°C.  We  might  also  permit  larger  standard  deviations 
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Y,  V, 

Food  Number  Initial  Initial 

and  name  tamp,  °C  count 

1 chicken  5 6 

2 Sauce  5 6 

3 Stock  and  Seasoning  5 6 

5 Chilled  vegetables  5 7 

7 Frozen  peas  -18  5 

9 Slurry  10  6 

12  Topping  5 6 


Table  1:  Initial  Attributes  of  Foods 


W, 

Weight,  lbs 

60 

2 

23 

8 

6 

2 

5 
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Figure  3:  Graphs  of  Y (Temperature  in  °C)  and  V (Natural  Logarithm  of  Bacteria  Count) 

versus  1 (time  step)  in  Foods  1,  4,  6 and  8 


in  some  of  the  random  numbers,  like  those  governing  the  storage  temperatures  and  the 
handling  errors.  In  this  way  we  could  see  how  large  the  bacteria  counts  become  as  we 
let  the  process  become  more  sloppy. 

It  is  perhaps  worth  emphasizing  that  the  two  monitor-systems  give  signals  when  the 
bacteria  count  is  high  enough,  regardless  of  how  close  the  food  is  to  being  served.  Food 
that  is  highly  contaminated  and  about  to  be  served  is  certainly  more  dangerous  than  an 
equally  contaminated  food  that  will  later  be  cooked  before  it  is  served.  However,  the 
two  monitoring  systems  ignore  that  distinction.  It  would  not  be  too  difficult  to  change 
the  monitor  systems  so  they  only  operated  on  foods  that  are  served  or  so  they  became 
increasingly  sensitive  as  the  process  went  on. 

It  is  evident  that  the  program  as  described  here  is  workable.  Undoubtedly 
improvements  or  extensions  of  it  are  possible,  but  it  is  necessary  to  try  it  first  on  a 
great  variety  of  processes  to  see  what  changes  need  to  be  made.  The  author  hopes  that 
it  will  be  convenient  and  useful  in  its  present  form. 
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Appendix  I:  Formulas  Usad  in  the  Simulation 


In  this  Section  we  give  more  complete  descriptions  of  the  calculation  of  bacterial 
change  during  storage  and  at  a handling  operation.  Also  we  describe  the  (wo  monitoring 
systems  more  thoroughly  than  before. 

A.  Bacterial  Growth  During  Storage  of  a Food 

We  assume  that  the  differential  equation  governing  the  non-random  growth  of 
organism  population  is 

dp/dt  = pa[Y(t)] 

where  p = density  of  organisms  (i.e.  number  of  organisms  per  gram),  t * time,  and  a[Y(t)] 
is  a function  that  depends  on  the  temperature,  Y,  which  itself  depends  on  t.  This  equation 
can  be  written 

dV/dt  = cr  [ Y (t)] 


if  we  define 


V = log8  p. 

Integrating  this,  we  obtain 

t 

V(t)  = V(tp)  + / a[Y(s)]ds, 


where  tp  is  any  time  earlier  than  t.  We  now  change  to  a discrete  approximation  in  which 
time  is  only  defined  at  t,,  tJ(  t3,  ...  where  for  all  i 

tj  - tM  = At. 

If  we  let  t * tj  and  tp  * tj -j,  abbreviate  V(tj)  by  Vj  and  V(tj_ -j ) by  and  use 

the  trapezoidal  approximation  for  the  integral,  we  obtain 
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Vj  = Vj.,  + 'h  At  n[ Y(tj)]  + a[Y(tj_1>] 


Then,  defining  for  all  i 

Aj  = At  Of [ Y (tj)] , 


we  get 


Vj  = VM  + % (Aj  + Aj.,,. 

We  add  to  this  the  random  change,  rj,  to  obtain  the  formula 
Vj  = Vj -)  + V4  (Aj  + Aj.,)  + rj, 

where  rj  are  uncorreiated  Gaussian  random  variables  with  mean  mr,  standard  deviation 


We  now  write  down  the  specific  formulas  defining  Aj.  The  temperature  of  a food 
is  found  from 


Yj  * Ys  + (Y,  - Ys)e“P'(i-'*) 

where  Y{  is  the  ambient  temperature  in  the  storage  area,  i is  the  present  time-step,  ij 
is  the  time-step  at  which  the  food  entered  the  storage  area,  Y,  is  the  temperature  of 
the  food  at  time  i]  and  P,  is  the  effective  thermal  diffusion  coefficient  of  the  food. 
Aj  is  expressed  in  terms  of  Yj  by  means  of 

Aj  = 0 if  Yj  < P, 

Aj  = P3  exp  -P4(P5-Yj)J  if  P,  < Y;  < Ps 

A;  - Pj-PftlYj-P,)3  if  Yj  / P, 

In  these  formulas  P2  is  the  temperature  below  which  (it  is  assumed)  there  is  negligible 
growth  or  death.  P3  is  the  temperature  at  maximum  growth,  and  P3  is  the  maximum 
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value  of  Aj,  i.e.  the  value  of  Aj  when  Y;  ■ P*.  P4  determines  how  rapidly  Aj  decreases 
from  P,  when  Yj  decreases  from  Pf.  P6  determines  how  quickly  Aj  decreases  from  P3 
as  Yj  increase*  above  P*.  When 

Yj  > P5  + iP3/P«)1/3, 

Aj  becomes  negative,  and  |Aj|  is  proportional  to  death  rate  rather  than  growth  rate. 

The  following  relationship  is  sometimes  useful  in  estimating  the  value  of  A.  Suppose 
that  a food  is  held  at  a constant  temperature  such  that  A > 0,  long  enough  so  that 
the  number  of  bacteria  is  just  double  its  initial  value,  i.e.  V has  increased  by  In  2 ■ .693. 
If  n time-steps  are  required  to  accomplish  this,  then  the  (constant)  value  of  A is 

A = (In  2 )/n  * .693/n 

and  nAt  is  the  generation  time,  usually  called  gj . Similarly,  if  A < 0 and  n is  the  number 
of  time  steps  needed  to  reduce  the  bacteria  count  to  half  its  initial  value,  then 

A = -.693/n 

and  nAt  is  the  half-life  of  the  organism  population. 

B.  Attributes  of  Quantities  Created  During  Handling 

Here  we  record  the  formulas  used  by  the  subroutine  SrLEND  in  deriving  the 
temperature,  Y,  weight,  W,  and  bacteria  count,  V,  of  new  foods  created  at  a split  or 
blend. 

At  a split  a single  food,  which  exists  and  has  attribute  values  Yq,  Wq,  Vq,  disappears, 
and  N new  foods  are  created  with  attribute  values  Y|<,  W|<,  V«  where  K * 1,  2 ...  N. 
We  define 

*K  * WK/WQ, 


i.e.  f|(  are  the  proportions  by  weight  of  the  new  foods,  and 
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f.  + fj  + ...  + fN  = 

These  proportions  are  specified  by  the  user.  Then  aside  from  handling  errors  the  attributes 
of  the  new  foods  are  taken  to  be 

YK  - Yo 

V'k  - v0  K - 1.  2.  ...  N 
WK  - *KW0 


Handling  errors  are  assumed  to  occur  in  V«  and  and  to  have  the  form 
VK  = VK  + lrVKl 
WK  = Wk  (1  lrWKl) 

where  ryj^  and  ryy are  Gaussian  — distributed  random  numbers  with  zero  mean  but 
different  variances.  The  nature  of  these  errors  follows  from  the  assumption  that  handling 
the  food  can  result  only  in  increased  counts  (i.e.  handling  contamination)  and  decreased 
weight  (spillage). 

At  a blend  N foods  with  attributes  Y^,  WK  and  VK,  K = 1,  2,  ...  N,  disappear, 
and  a single  food  with  attribute  values  Yq,  and  Wq  and  Vq  is  created.  Disregarding 
handling  errors  temporarily,  we  calculate  Wq,  Yq  and  Vq  as  follows.  The  weight  of 
the  new  food  is  the  sum  of  the  weights  of  the  foods  being  blended,  i.e. 

Wq  = W,  + W3  + ...  + WN  = ZWK 

The  temperature  of  the  new  food  is  the  weighted  average  of  the  temperatures  of  the 
foods  entering  the  blend,  i.e. 

Yo  = W'Y>  + WjY*  + •••  < wnyn  = 2Wkyk 

W,  + Wj  + ...  + WN  Wq 
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To  calculate  Vq,  we  recall  that  the  density  of  bacteria  in  the  K-th  food  is  p«, 


PK  * «VK. 

and  the  total  number  of  bacteria  in  the  K-th  food,  is  the  density  multiplied  by  the 
weight, 

-K  ' *KPK  * 

The  total  number  of  bacteria  in  the  new  food,  vq,  is  the  sum  of  the  numbers  in  the 
foods  entering  the  blend, 

vQ  * 2>K. 

Finally  Vq  is  the  natural  logarithm  of  the  density  of  organisms  in  the  new  food,  i.e. 
Vo  - loge  ( vq/VJq  } = log,,  . 

L 2wkJ 

Handling  errors  in  a blend  are  assumed  to  have  the  same  form  as  in  a split,  i.e. 

vo  * vq  + |rVQ| 

wO  ■ wo 

where  r^Q  and  r^Q  are  Gaussian  - distributed  random  numbers  with  zero  mean  but 
(in  general)  different  variances. 

C.  The  Two  Monitor  Systems 

The  record-keeping  monitor  system  operates  as  follows:  each  food  is  assigned  an 
estimated  bacteria  count,  RI(J),  when  it  first  comes  into  the  system.  Each  food  leaving 
the  system  before  the  end  of  simulation  is  also  assigned  a final  count  RO(J),  in  a manner 
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to  be  described  shortly.  For  a food  that  is  in  the  system  at  the  start  of  simulation, 
or  is  brought  in  from  outside  during  simulation,  the  value  of  R I (J)  is  assigned  by  the 
user,  based  on  previous  experience  with  that  type  of  food.  Foods  created  out  of  other 
foods  during  simulation  get  their  values  i><  R I { J ) from  the  RO(J)  values  of  their  progenitor 
foods  in  the  following  ways:  if  the  JO-th  food  comes  from  a split  in  the  J-th  food 

R I (JO)  = RO(J) 

If  the  JU-th  food  comes  from  a blend  of  foods  numbered  JR  I,  JR  I,  ....  JRN,  then 

RI(JU)  = 1 (RO(JRI)  + ROJR2)  + ...  + RO(JRN)]. 

N 

We  now  describe  how  RO(J)  is  determined.  Three  parameters,  Ql,  Q2,  and  Q3, 
are  supplied  by  the  user,  and  a temperature,  YCOMP,  is  taken  to  be  the  highest  of  all 
the  ambient  temperatures  to  which  the  J-th  food  is  exposed.  Then 

R 0( J)  = R I ( J)  if  YCOMP  < Ql 

= R I ( J)  + Q3  (YCOMP-QIMIT-II) 
if  Ql  < YCOMP  < Q2 

= O if  YCOMP  > Q2. 

In  these  formulas  Ql  is  a temperature  below  which  no  growth  or  death  occurs,  Q2  is 
a temperature  above  which  the  organism  is  almost  completely  killed  off  II  and  IT  are 
the  time-steps  at  which  the  food  enters  and  leaves  the  system  respectively,  and  Q3  is 
a coefficient  specifying  the  estimated  bacterial  growth  rate. 

In  the  program  the  value  of  RO{J)  for  each  food  is  compared  with  the  number  TESTI, 
supplied  by  the  user.  TESTI  is  the  level  of  contamination  that  the  monitor  system  considers 
to  be  "bad".  If  RO(J)  exceeds  TESTI,  this  constitutes  a signal  that  the  monitor  system 
has  detected  an  unsafe  situation  in  that  food.  A message  to  this  effect,  "MON R EC 
SIGNAL",  will  appear  in  the  printed  output  for  the  food.  Also,  if  several  repetitions 
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of  the  simulation  are  made,  the  program  prints  out  NS,  the  number  of  repetitions  in 
which  this  monitor  system  detected  an  unsafe  situation. 

Finally  a number  called  TRT,  supplied  by  the  user,  is  the  bacteria  count  which  he 
takes  as  the  actual  count  at  which  a "bad"  (dangerous)  condition  occurs.  This  need 
not  be  the  same  as  TESTI,  the  number  which  the  monitor  system  is  using  for  that  purpose. 
The  program  calculates  the  largest  bacteria  count  occurring  during  the  simulation  and 
compares  it  with  TRT,  to  determine  whether  a dangerous  situation  actually  does  occur 
during  the  simulation.  If  several  repetitions  of  the  simulation  are  done,  the  program  prints 
out  NT,  the  true  number  of  trials  in  which  a dangerous  situation  occurres.  To  help 
in  assessing  the  effectiveness  of  the  monitor  system,  two  other  numbers,  Nil  and  N00, 
are  also  printed  out.  Nil  is  the  number  of  trials  in  which  the  monitor  system  correctly 
stated  that  there  was  a dangerous  situation,  and  N00  is  the  number  of  trials  in  which 
the  monitor  system  correctly  stated  that  there  was  no  dangerous  situation. 

The  domestic  monitor  system  operates  quite  differently.  It  is  assumed  that  at  each 
food-handling,  i.e.  split  or  blend,  a cook  checks  the  bacteria  count  of  the  food  coming 
into  the  split  or  emerging  from  the  blend.  Conceptually,  we  imagine  that  the  cook  sniffs 
the  food  and  detects  a signal  if  the  food  smells  "bad".  In  the  program  this  is  done 
by  the  subroutine  SNIFF,  which  compares  the  bacteria  count,  VT,  with  a random  number, 
C,  that  represents  the  cook's  threshold  for  the  detection  of  bad  smells.  If  VT  exceeds 
C,  a detection  is  made.  The  threshold,  C,  is  obtained  from 

C = AV  + SD'XX 

where,  XX,  is  a standardized  Gaussian  random  number.  Thus,  AV  represents  the  average 
threshold  value  and  SO  the  standard  deviation  in  the  threshold  value. 

There  is  a further  aspect  to  the  domestic  monitor  system.  It  is  assumed  that  at 
each  food  handling,  the  cook  goes  to  all  the  storage  units  containing  one  or  more  of 
the  foods  involved  in  the  split  or  blend.  He  then  checks  all  the  foods  resident  in  those 
storage  areas  at  that  time.  In  the  program  this  is  done  by  first  determining  the  storage 
areas  involved,  then  finding  what  foods  are  in  these  storage  areas  and  finally  applying 
the  subroutine  SNIFF  to  these  foods. 
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It  is  reasonable  to  expect  that  the  cook's  threshold  for  detection  would  be  higher 
in  the  situation  just  described  than  in  the  previous  circumstances,  where  he  is  checking 
food  that  he  is  actually  working  on.  Accordingly  the  values  of  A V and  SD  are  permitted 
to  be  different  in  the  two  different  applications  of  SNIFF.  AVI  and  SD1  are  the 
parameters  of  the  threshold  random  number  for  foods  being  directly  handled,  and  AV2 
and  SD2  are  the  corresponding  parameters  when  the  foods  are  being  only  indirectly 
checked. 

When  there  are  several  repetitions  of  the  simulation,  the  program  prints  out  quantities 
NSD,  NDOO  and  ND11,  analogous  to  the  NS,  N00  and  Nil  printed  out  for  the 
record-keeping  monitor  system. 
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Appendix  II:  Input  Data 


In  discussing  the  input  data  format  we  shall  always  describe  it  as  if  it  were  stored 
on  cards,  because  this  is  the  storage  medium  used  in  entering  data  for  the  first  time. 

The  general  arrangement  of  the  input  data  is  as  follows.  First,  there  are  four  cards 
giving  general  information  about  the  simulation,  such  as  NJ,  the  number  of  separate  foods, 
NMS,  the  number  of  splits,  NSTOR  the  number  of  storage  units,  as  well  as  information 
about  the  two  monitor  systems.  Next  are  the  cards  describing  the  properties  of  the 
individual  foods.  Each  food  requires  at  least  (usually  exactly)  six  cards  to  define  it. 
Following  these  is  a single  card  specifying  the  temperatures  of  the  storage  units.  Then 
come  the  cards  supplying  information  about  the  splits,  three  cards  per  split,  and  finally 
the  cards  giving  the  data  for  the  blends,  again  three  cards  for  each  blend. 

There  are  only  two  essentially  different  formats  for  input  cards.  Most  of  the  cards 
are  in  free-field  format,  i.e.  they  consist  of  numbers  (some  with  decimal  points,  some 
without)  separated  by  commas,  and  without  any  blanks.  All  the  remaining  cards  are 
of  "Legend"  type,  that  is,  they  are  messages,  made  up  of  ordinary  words,  used  to  annotate 
the  printed  output  of  the  simulations. 

It  is  necessary  to  be  careful  about  the  order  in  which  the  individual  foods  are  entered. 
The  program  assigns  J=1  as  the  number  of  the  first  food  entered,  J=2  that  of  second 
food  and  so  on.  The  main  body  of  the  program  consists  of  calculations  done  in  the 
following  order: 

(a)  BACGRO  for  J=1  (i.e.  bacterial  growth  of  food  number  1 during 
storage) 

(b)  SPLEND  for  J=1  (i.e.  split  or  blend  whose  input  is  J=1,  if  it  can  be 
carried  out) 

(c)  - BACGRO  for  J=2 

(d)  SPLEND  for  J=2 
etc. 
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Because  of  this  arrangement  the  foods  must  be  entered  in  such  an  order  that  foods  created 
at  splits  or  blends  are  entered  after  their  progenitor  foods.  For  example,  if  cake  flour 
is  created  by  a split  in  flour,  and  the  data  for  cake  flour  were  inadvertently  entered  before 
flour,  then  the  program  would  assign  J=1  to  cake  flour  and  J=2  to  flour.  The  subroutine 
SPLEND  would  then  wrongly  derive  the  initial  bacterial  count  and  temperature  for  flour 
from  those  of  cake  flour. 

This  requirement,  that  foods  created  at  a split  or  blend  be  entered  after  their 
progenitor  foods,  is  not  a severe  one  and  tends  to  agree  with  what  a user  would  ordinarily 
do  in  any  case.  The  main  advantage  of  preparing  a diagram  like  Figure  2 is  that  one 
can  easily  check  that  this  condition  is  satisfied  by  making  sure,  that  the  numbers  of  the 
foods  on  the  right  of  (emerging  from)  each  split  or  blend  are  higher  than  those  on  the 
left  (entering).  When  the  diagram  is  such  that  this  condition  is  satisfied  at  all  the  splits 
and  blends,  then  entering  the  foods  from  top  to  bottom  of  the  diagram  will  give  a correct 
order. 

Incidentally,  an  alert  reader  will  notice  that  in  Figure  2 all  the  foods  entering  a 
blend  are  numbered  consecuti  'ely.  This  condition  makes  the  diagram  clear  and  easy  to 
read,  but  the  program  does  not  require  it.  I.e.,  the  program  will  function  correctly  even 
though  the  foods  entering  a blend  are  not  numbered  consecutively,  provided  the  previous 
condition  is  met. 

We  now  describe  the  information  on  each  card,  and  give  as  an  example  the 
corresponding  information  for  the  chicken  pot  pie  simulation  of  Section  4.  The  cards 
are  grouped  according  to  the  general  arrangement  stated  near  the  beginning  of  this 
Appendix.  Following  this  detailed  description  of  the  input  is  a printed  listing  of  all  the 
input  data  for  the  simulation  of  chicken  pot  pie. 
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Group  A:  General  Information 


Card  1 : Legend  Card.  Any  desired  message  can  be  punched  in  columns  2-72  and 

will  be  printed  at  the  top  of  the  first  output  page. 

Example:  CHICKEN  POT  PIE  SIMULATION 

is  punched  in  columns  2-27 

Card  2:  NJ,  NMS,  NMB,  X(1),  IMAX,  SDV,  NREP,  NRNO,  IPTRD,  NSTOR 

Example:  13,  1,  4,  347.,  100,  .05,  5,  130,  0,  4 
is  punched  in  columns  1-29 

Here 

NJ  = number  of  separate  foods  * 13 
NMS  = number  of  splits  *1 
NMB  = number  of  blends  = 4 

X(1)  = the  "seed"  for  random-number  generation  = 347. 

IMAX  = an  integer  larger  than  the  total  number  of  time-steps  in  the 
simulation  * 100 

SDV  * standard  deviation  of  the  random  error  in  the  bacterial  growth 
equation  * .05 

NREP  = number  of  repetitions  of  the  simulation  « 5 
NRNO  = number  of  random  numbers  generated  at  each  repetition  = 130 
IPTRD  = controller  for  printing.  If  IPTRD  * 1,  printing  of  Y,  V and  YA 
is  suppressed  except  at  the  initial  and  final  times  for  each  food. 
Choosing  IPTRD=0  in  the  example  permits  printing  of  all  values 
of  Y,  V and  YA. 

NSTOR=  number  of  storage  units  = 4 

Card  3:  AVI,  SD1,  AV2,  SD2 

Example:  12.,  1.,  13.,  1 

is  punched  in  columns  1—13 
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See  the  discussion  of  domestic  monitor  system.  Appendix  I.C. 


AVI  = average  threshold  count  for  directly-handled  foods  * 12. 

SD1  = standard  deviation  in  these  threshold  counts  = 1. 

AV2  = average  threshold  count  for  foods  not  directly  handled  * 13. 
SD2  = standard  deviation  in  these  threshold  counts  * 1. 

Card  4:  Q1,  Q2,  Q3,  TESTI,  TRT 

Example:  10.,  50.,  .01,  12.,  12. 

is  punched  in  columns  1-19 

See  the  discussion  of  the  record-keeping  monitor  system.  Appendix  I.C. 

Q1  = lower  end  of  temperature  growth  range  = 10. 

Q2  = upper  end  of  temperature  growth  range  = 50. 

Q3  = bacterial  growth  rate  = .01 

TESTI  = threshold  of  dangerous  contamination  used  by  monitor  = 12. 
TRT  = true  threshold  of  dangerous  contamination  = 12. 

NOTE  1:The  simulation  assumes  that  NMS  > 1 and  NMB  > 1,  see  Group  A,  Card  2, 
i.e.  there  is  a least  one  split  and  one  blend.  These  must  be  introduced  if  they  do  not 
occur  naturally.  For  example,  if  one  wants  to  simply  find  the  growth  of  bacteria  in 
chicken,  then  one  must  introduce  a one-for-one  split  and  one-for-one  blend,  giving  the 
chicken  three  different  food  numbers  in  the  simulation.  This  means  that  we  must  always 
have  NJ  > 3. 

Group  B:  nformation  specific  to  a Food. 

For  each  food  the  following  set  of  cards  must  be  prepared.  The  examples  given  are 
for  tf  first  food  in  the  chicken  pot  pie  simulation. 

Caro  1:  Legend  Card.  Any  desired  description  or  message  can  be  punched  in 

columns  2-72  and  will  be  printed  along  with  the  value  of  J at  several 
points  in  the  output. 
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Example:  CHICKEN 

is  punched  in  colums  9—15 

Card  2:  IC(J),  W(J),  Y(J,1),  V(J,1),  ISER(J),  MST(J) 

Example:  1,  60.,  5.,  6.,  100,  3 

is  punched  in  columns  1-17 

IC(J)  = the  time  step  at  which  food  is  first  in  the  system  * 1 

W(J)  = the  intended  weight  of  the  food  * 60, 

Y(J,1)  * the  temperature  at  the  initial  time  step  * 5. 

V(J,1)  = the  initial  bacteria  count  * 6. 

ISER(J)  * is  an  integer  which  is  the  time-step  at  which  the  food  is  served. 

If  the  food  is  not  served  (as  in  this  example)  ISER(J)  must 
be  at  least  as  large  as  (MAX,  ISER(J)  = 100 
MST(J)  = the  storage  unit  w^iere  this  food  is  stored  - 3 

Card  3:  ICD(J),  WD(J),  YD(J),  VD(J),  NATC(J),  CONT(J) 

Example:  0,  0.,  0.,  0.,  1,0 

is  punched  in  columns  1—15 

ICD(J)  * standard  deviation  of  random  error  in  initial  time  * 0. 

WD(J)  * standard  deviation  of  random  error  in  weight  * 0. 

YD(J)  * standard  deviation  of  random  error  in  initial  temperature  = 0 
VD(J)  = standard  deviation  of  random  error  in  initial  bacteria  count  = 0. 

NATC(J)  = the  number  of  ambient  temperature  changes  in  the  food  * 1 

CONT(J)  = the  mean  of  the  random  increase  in  bacteria  count  = 0. 

Card  4:  P(J,1),  P(J,2),  P(J,3),  P(J,4),  P(J,5),  P(J,6) 

Example:  .17,  0.,  .15,  .0047,  40.,  .000443 

is  punched  in  columns  1—28 
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For  the  J-th  food  these  are  the  constants  referred  to  as  Pi,  Pj, ...  P6  in  Appendix  I.A. 
They  describe  approximately  the  behavior  of  E.  coli  at  five-minute  time  steps. 


P(J,1) 

P(J,2) 

P(J,3) 

P(J,4) 

P(J,5) 

P(J,6) 

= thermal  diffusion  coefficient  = .17 
= temperature  below  which  no  growth  or  death  occurs  * 0. 
= maximum  value  of  growth  rate  = .15 
= coefficient  of  growth  rate  change  below  maximum  growth 
temperature  = .0047 

= temperature  of  maximum  growth  rate  = 40. 

= coefficient  of  growth  rate  change  above  maximum  growth 
temperature  = .000443 

Card  5: 

R 1 (J) 

Example: 

6 

is  punched 

in  columns  1-2 

RHJ) 

= initial  bacteria  count  assumed  by  record-keeping  monitor 

= 6. 

Card  6: 

IAC(J,1 ),  IACD(J,1),  YAC(J,1) 

Example: 

4,  0,  100 

is  punched  in  columns  1-8 


IAC(J,1)  = time  step  at  which  first  ambient  temperature  change  occurs 
= 4 

IACD(J,1)  = stjndard  deviation  of  the  random  error  in  the  time  step  at 
which  first  ambient  temperature  change  occurs  = 0 
YAC(J,1)  = the  ambient  temperature  after  its  first  change  = 100. 

Note  1 : If  more  than  one  ambient  temperature  change  occurs,  one  additional  card,  exactly 
like  card  6,  should  be  made  up  for  each  additional  temperature  change  and  added  following 
card  6.  The  cards  describing  ambient  temperature  changes  should  be  in  order  of  increasing 
values  of  IAC. 


27 


Not*  2:  If  no  ambient  temperature  change  occurs,  card  6 should  contain  the  entries 
0,  0,  0.  in  columns  1-6.  Card  6 is  then  the  final  card  for  that  food. 

Note  3:  For  foods  not  created  during  the  simulation  the  program  uses  all  the  information 

on  these  input  cards.  For  foods  created  during  the  simulation  the  quantities  W(J),  Y(J,1) 
and  V(J,1)  on  card  2 of  this  group  are  not  used.  Instead  the  program  obtains  these 
quantities  from  the  attributes  of  the  progenitor  foods  by  using  subroutine  SPLEND.  Hence 
any  reasonable  values  can  be  entered  for  these  quantities,  in  particular  0.,  0.,  0. 

Group  C:  Information  about  Storage  Units 

This  group  contains  a sing'e  card  with  2*NSTOR  numbers,  giving  the  intended  ambient 
temperature  and  standard  deviation  of  the  random  error  in  that  temperature  for  each 
of  the  NSTOR  storage  units. 

Card  1:  YST(1),  YSTD{1),  YST(2),  YSTDI2) YST(NSTOR),  YSTD(NSTOR) 

Example:  -18.,  0.,  5.,  0.,  20.,  5.,  100.,  10 

is  punched  in  columns  1—29 

YST(1)  * ambient  temperature  in  storage  unit  1 * -18. 

YSTD(1)  * standard  deviation  of  random  error  in  that  temperatir  = 0. 
YST  (2)  = ambient  temperature  in  storage  unit  3*5. 

YSTD(2)  = standard  deviation  of  random  error  in  that  temperature  = 0. 
etc. 

Group  D:  Information  about  Splits 

The  split  number,  MS,  is  assigned  to  each  split  in  order  as  its  data  is  read  in.  Exactly 
three  cards  define  the  data  of  each  split.  See  Appendix  I.B. 

Card  1:  NJO(MS) 

Example:  2 

is  punched  in  column  1 

NJO(MS)  = the  number  of  different  foods  created  at  split  number  MS  * 2 
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Card  2:  IS(MS),  JSI(MS),  ISD(MS).  JS0(MS,1),  PR(MS,1),  JS0(MS,2),  PR(MS,2), ... 


Example:  13,  1,  0,  4,  .3,  10,  .6 

IS(MS)  = time-step  at  which  split  number  MS  occurs  = 13 

JSI(MS)  = number  of  the  food  being  split  = 1 

ISD{MS)  = standard  deviation  of  the  random  error  in  the  time  of  the 

split  = 0 

JSO(MS,1)  = number  of  the  first  food  created  at  the  split  = 4 

PR(MS,1)  = proportion  (by  weight)  of  the  first  food  created  at  the  split 

= .3 

JSO(MS,2)  = number  of  the  second  food  created  at  the  split  = 10 
PR(MS,2)  = proportions  of  the  second  food  created  = .6 

Card  3:  HEVS(MS,1),  HEWS(MS,1),  HEVS(MS,2),  HEWS(MS,2)  ... 

Example:  .3,  0.0,  .5,  .05 

HEVS(MS,1)  = standard  deviation  of  the  random  handling  contamination  in 
the  first  food  created  at  the  split  = .3 

HEWS(MS,1)  = standard  deviation  of  the  random  spillage  weight  loss  in  the 
first  food  created  at  the  split  = 0.0 

HEVS(MS,2)  = standard  deviation  of  contamination  in  the  second  food 
created  at  the  split  = .5 

HEWS(MS,2)  = standard  deviation  of  spillage  loss  in  the  second  food  created 
at  the  split  = .05 

NOTE:  The  time  of  the  split,  IS(MS),  must  always  agree  with  the  initial  times,  IC(J), 

of  all  the  foods  created  at  the  split.  Also  ISD(MS)  must  agree  with  ICD(J)  for  all  the 
foods  created  at  the  split,  and  the  latter  must  agree  with  each  other. 
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Group  E:  Information  about  Blands 


The  blend  number  is  assigned  to  each  blend  in  order  as  its  data  is  read  in.  Exactly 
three  cards  define  the  data  of  each  blend.  Sec  Appendix  I.B. 

Card  t:  NJI(MB) 

Example:  4 

is  punched  in  column  1 

NJI (MB)  * the  number  of  different  foods  blended  together  at  blend 
MB  - 4 


Card  2:  IB(MB),  JBO(MB),  IBD(MB),  JBI(MB,1),  JBI(MB,2),  JBI(MB,3),  JBI(MB,4) 

Example:  17,  6,  0,  2,  3,  4,  5 


IB(MB) 

JBO(MB) 

IBD(MB) 

JBI(MB,1) 

JBI(MB,2) 

etc. 


time  step  at  which  blend  number  MB  occurs  * 17 
the  number  of  the  food  created  at  the  blend  * 6 
standard  deviation  of  the  random  error  in  the  time  of  the 
blend  * 0 

number  of  the  first  food  blended  * 2 
number  of  the  second  food  blended  * 3 


Card  3:  HEVB(MB),  HEWB(MB) 

Example:  .4,  .02 

HEVB(MB)  - standard  deviation  of  the  random  handling  contamination  in 
the  food  created  at  blend  MB  - .4 

HEWB(MB)  * standard  deviation  of  the  random  spillage  weight  loss  in  the 
food  created  at  the  blend  * .02 
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NOTE:  The  time  of  the  blend,  IB(MB),  must  always  agree  with  the  initial  time,  IC(J), 
of  the  food  created  at  the  blend.  Also  IBD(MB),  its  standard  deviation,  must  agree  with 
ICD(J)  for  the  food  created  at  the  blend. 

This  concludes  the  description  of  the  input.  Below  is  a complete  listing  of  the  input 
data  for  the  chicken  pot  pie  simulation.  Each  line  of  printing  represents  the  data  punched 
on  one  input  card. 
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1 chicken  pot  pie  simulation 

! 13»l»4»347.»itO*,Q»*0*WQ»Q*4 
6.09*0.* 7*!» 0.  — ^ 

io.»9o.*.u*s.*7,i9 
CHICKEN  . \ 

1*60. *5. *64*100*1 
0# 0. » 0* » 0* * 1*0.  \ 

, . 17*0, * .15* .0047*40. » *000443 
6 • 

4*0*100. 

SAUCE  (ROUX) 

1*2. *5«*6. *100*2 
0 *0 . * 0. »0, * 1*0. 

.50*0.*. 15*. 0047*40. *.000443 

6. 

13*1*100.  1 

STOCK  ANO  SEASONING 
1 *23.  *5. *6. *100*2 

0f0.*0.*0.*0*0. 

•50*Q.* .15* .0047*40. *.000443 
6 . 

0*0*0. 

extra  stock 

13*0. *0. *6. * 100  * 3 
1*0. *0. *0. *0*0. 

• 50*0.  *.15* .0047*40. *.000443 

6. 

0*0*0. 

CHILLEO  VEGETABLES 
1*G.*S.»7. *100*2 
0*0.* 0**0. *0*0. 

• 50*0.*. II*  .0047(40*  *.000443 

6. 

0*0*0. 

MIXTURE  ALPHA 
17* 0.»0.»0. *100*4 
! 1*0. »0.*0. *0*0. 

•50*0.* .11* .0047*40. * .000443 

6t 

r 0*0*0. 

FROZEN  PEAS 
1 *G. * *10. *3. *100* 1 

o»o.*o.*o.*o»o. 

.50*0. * .19*. 0047*40. * .000443 

6. 

0*0*0. 

MIXTURE  BETA 
25*0.»0.*0. *100*4 
1»0.»0.*0.*0*0. 

•40*0. » .11* .0047*40. *.000443 
0* 

1 0*0*0. 

I slurry 

25 *2. * 10. *0. *100 *3 
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0*0«*0.*0.»1*0. 

•40*0. » .11 # .0047*60. * .000443 

6. 

aT»o#ioo. 

OI-«OrCO»  COOKED  CHICKEN 

1 * 0. » 0 • *0* * l • •# 

.17*0, *.15*. 0047*48 .*.000443 

6. 

27*0*20. 

MIXTUKE  0AMMA 
31*0. *0.*0. *108*4 
1*0.*8**0.*#*8. 

• 30* 0, » .15*  .0047*40. » • 000443 

6. 

0*0*0* 

topping 

35*5. *5. *6. *100*4 
0 » 0 . * 0* *0. *0*0 • 

.17* 0.*. 15 *.0047 *40. *.000443 
6. 

0*0*0. 

CHICKEN  POT  PIE 
39  *0.*0.*0.i 97  * 3 
2*0*  *0.»0. *2*0. 

.20*0. *.15* .0047*40* *.000443 
6 • 

45*2*-18» 

91*2*175. 

«18.»0.*5«»0.»20.*5.*100.*10. 

2 

13* 1*0*4*. 3*10*. 6 
.3*0.0* .5* .05 
4 

17*6*0*2*3*4*5 
•4* .02 
2 

25*8*0*6*7 
j 3 » .01 
3 

31*11*0*8*9*10 
. 3 * . 05 
2 

39*13*0*11*12 
.4* .04 
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Appendix  III:  Output 


In  this  Appendix  we  discuss  the  printed  output  from  the  simulation  of  chicken  pot 
pie  described  in  Section  4.  The  general  form  of  the  output  is  sketched  in  the  following 
paragraphs,  and  the  actual  print-out  from  a small  simulation  is  shown.  The  Appendix 
concludes  with  some  detailed  comments  about  this  output. 

The  output  can  be  divided  into  two  main  sections.  The  first  contains  a listing  of 
all  the  information  that  the  user  furnishes  about  the  simulation,  i.e.  a methodical  record 
of  the  input  data.  Within  this  section  there  are  a number  of  subsections,  conforming 
with  the  grouping  of  input  data  in  Appendix  II.  The  first  subsection  contains  the 
descriptive  title  of  the  run,  followed  by  four  lines  of  general  information.  The  data 
pertaining  to  the  individual  foods  is  listed  next,  consisting  of  an  annotative  heading  plus 
(usually)  five  lines  of  data  about  each  food.  The  third  subsection  presents  the  data  on 
the  storage  units,  and  the  last  subsection,  headed  "Intended  Schedule  of  Splits  and  Blends", 
summarizes  the  information  about  splits  and  blends. 

The  second  main  section  describes  the  actual  (as  contrasted  with  intended)  events 
of  the  simulation  and  its  results.  It  begins  with  the  table  headed  "Final  Schedule",  followed 
by  a summary  of  information  about  each  food  after  the  introduction  of  random  errors. 
The  next  subsection  contains  the  calculated  time-dependence  of  Y,  food  temperature,  V, 
logarithmic  bacteria  count.  A,  growth  rate  and  YA,  storage  temperature,  for  each  food 
during  its  storage,  followed  by  what  happens  to  the  food  at  the  end  of  storage.  This 
subsection  contains  the  principal  results  of  the  simulation.  Each  repetition  of  the 
simulation  leads  to  another  section  of  output  starting  with  the  "Final  Schedule"  table. 
After  the  output  from  all  repetitions  is  finished,  there  is  a final,  two-line,  subsection 
summarizing  the  effectiveness  of  the  two  monitor  systems. 

At  this  point  it  is  useful  for  the  reader  to  examine  the  output  for  the  sa  nple  run 
of  the  chicken  pot  pie  simulation  whose  input  was  listed  at  the  end  of  Appendix  2.  The 
reader  will  notice  in  the  line  after  the  heading  that  NREP=5,  i.e.  this  run  included  5 
repetitions  of  the  simulation.  In  this  report  only  the  printout  from  the  first  and  fourth 
repetitions  are  shown  in  order  to  save  space. 
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SAQD.PE  MIMA6E.CHIP0TPIE 
CHICKEN  POT  PIE  SIMULATION 

NJ=13  NMS=  1 NMB=  4 X(i)s  347,  IMAX=  100  SQV=  ,050000 
NRNO=130  IPTRU=0  NSTORs  4 

AVIs  6*050  SD1=  .000  AV2=  7*200  S02=  .000 

Ql=  10.0000  Q2s  50.0000  Q3=  .0100  TEST1=  5.0000  TRT= 

J=  1 CHICKEN 

IC=  1 M=  60.0000  YI=  5.0000  Vis  6*0000  ISERslOO  STOR. 
ICDs  0 MDs  .0000  YD=  .0000  VD=  .0000  NATC=  1 CONT= 
P(J»K)  ,17000  .00000  .15000  .00470  40.00000 

Ks  1 IACs  4 IAC0=  0 YAC=  100.0000 
RIs  6.0000 


J=  2 SAUCE  (ROUX) 

IC=  1 Ms  2.0000  YI=  5.0000  Vis  6.0000  ISERslOO  STOR. 

ICDs  0 MDs  .0000  YDs  .0000  VDs  .0000  NATC=  1 CONT= 
P ( J , K ) .50000  .00000  .15000  .00470  40.00000 

K-  1 IAC=  13  1 ACDs  1 YAC=  100.0000 
RI=  6.0000 


J=  3 STOCK  AND  SEASONING 

ICs  1 M=  23.0000  Y I - 5.0000  Vis  6.0000  ISERslOO  STOR. 
ICDs  0 MDs  .0000  YD=  .0000  VD=  .0000  NATC=  0 CONT= 
P ( J r K ) .50000  .00000  .15000  .00470  40.00000 

Ks  1 IAC=  0 I ACDs  0 YAC=  .0000 
RI=  6.0000 


Js  4 

EXTRA 

STOCK 

ICs  13 

Ws 

.0000 

Yls  .0000  Vis 

6.0000 

ISERslOO 

STOR. 

ICDs  1 

MDs 

.0000 

YDs  .0000  VDs 

.0000 

NATC=  0 

CONTs 

P(J,K) 

.50000 

.00000  .15000  .00470  40 

.00000 

Ks  1 I ACs 

0 lACDs 

0 YACs  ,0000 

J=  5 

chilled  vegetables 

ICs  1 

Ms 

8.0000 

Yls  5.0000  Vis 

7.0000 

ISERslOO 

STOR. 

ICOs  0 

MDs 

.0000 

YDs  .0000  VDs 

.0000 

NATCs  0 

CONTs 

P(J.K) 

.50000 

.00000  .15000  .00470  40 

.00000 

Ks  1 I ACs 

0 lACQs 

0 YACs  .0000 

RIs 

6.0000 

Js  6 MIXTURE  ALPHA 

IC=  17  M=  .0000  YI=  .0000  Vis  .0000  ISERslOO  STOR. 

ICO-  1 WD=  .0000  YD=  .0000  VDs  .0000  NATOs  0 CONT= 

P ( J,K)  .50000  .00000  .15000  .00470  40.00000 

K=  1 IAC=  0 IACD=  0 YAC-  .0000 

Js  7 FROZEN  PEAS 

IC=  1 Ms  6.0000  YIs-18.0000  Vis  5.0000  ISERslOO  STOR. 

ICDs  0 MDs  .0000  YDs  .0000  VDs  .0000  NATC=  0 CONTs 

P ( J,K)  .50000  .00000  .15000  .00470  40.00000 

Ks  1 I ACs  0 XACDs  0 YACs  ,0000 
RIs  6.0000 


NREP=  5 


7.1500 


N0*=  3 
.0000 
•00044 


NO.=  2 
,0000 
.00044 


NO.=  2 
.0000 
.00044 


NO  • = 3 
.0000 
.00044 


NO.s  2 
.0000 
•00044 


N0.=  4 
.0000 
• 00044 


NO.=  X 
,0000 
, 00044 
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4z  G . ..  MIXTURE  fiEIA 

XC=  £5  Is  .0000  YI=  .0000  Vis  *0000  ISERslOO  STOR.  NO.=  4 

ICOl  I Mi  *0000  YDS  .0000  VDs  .0000  NATCS  0 CONTs  .0000 

P(J,K)  .40000  .00000  .15000  .00470  40.00000  .00044 

Ks  1 Ufii  0 IACDs  0 YACs  .0000 

J=  9 SLURRY 

IC=  25  Ws  2.0000  YU  10.0000  Vis  6.0000  ISERslOO  STOR.  NO.=  3 

ICOs  0 WDs  .0000  YD=  .0000  VDs  .0000  NATOs  1 CONTs  .0000 

P(J,K)  ,40000  .00000  .15000  .00470  40.00000  .00044 

Ks  1 IACs  27  IACD=  0 YACs  100.0000 

JslO  DE-BONED,  COOKED  CHICKEN 

IC=  13  Ws  .0000  YU  .0000  Vis  .0000  ISERslOO  STOR.  NO.s  2 

ICOs  1 WDs  ,0000  YDs  .0000  VDs  .0000  NATCs  1 CONTs  .0000 

P(J.K)  .17000  .00000  .15000  .00470  40.00000  .00044 

Ks  1 IACs  27  IACDs  0 YACs  20,0000 

jsll  MIXTURE  GAMMA 

ICs  31  Ws  .0000  YIs  .0000  Vis  .0000  ISERslOO  STOR.  NO.=  4 

ICOs  1 WDs  .0000  YDs  .0000  VDs  .0000  NATCs  0 CONTs  ,0000 

PU.K)  .30000  .00000  .15000  .00470  40.00000  .00044 

Ks  l IACs  0 IACDs  0 YACs  .0000 

Jsl2  TOPPING 

ICs  35  Ws  5.0000  YIs  5.0000  V Is  6.0000  ISERslOO  STOR.  NO.=  4 

ICDs  0 WDs  .0000  YDs  .0000  VDs  .0000  NATC=  0 CONT=  .0000 

P(J,K)  .17000  .00000  .15000  .00470  40.00000  .00044 

Ks  i IACs  0 IACDs  0 YACs  .0000 

Jsl3  CHICKEN  POT  PIE 

ICs  39  Ws  .0000  YU  .0000  V Is  .0000  ISERs  97  STOR.  NO.=  3 

ICDs  2 WDs  .0000  YDs  ,0000  VDs  .0000  NATCs  2 CONTs  .0000 

P ( J.K)  .20000  .00000  .15000  .00470  40.00000  .O0O44 

Ks  1 IACs  45  IACDs  2 YACs  -18.0000 
Ks  2 IACs  91  IACDs  2 YACs  175.0000 

STORAGE  UNIT  1 IS  AT  TEMP  -19.000  WITH  STD  DEV  .0000 

STORAGE  UNIT  2 IS  AT  TEMP  5,000  WITH  STD  DEV  .0000 

STORAGE  UNIT  3 IS  AT  TEMP  20.000  WITH  STD  DEV  5.0000 

STORAGE  UNIT  4 IS  AT  TEMP  100.000  WITH  STD  DEV  10.0000 

INTENDED  SCHEDULE  OF  SPLITS  AND  BLENDS 

SPLIT  NUMBER  1 IN  FOOD  NUMBER  1 OCCURS  AT  TIME  13  WITH  STD  DEV  0 
THE  NUMBERS  AND  PROPORTIONS  QF  THE  RESULTING  FOODS  ARE 
FOOD  NUMBER  4 PROPORTION  .30000 

FQOO  NUMBER  10  PROPORTION  .60000  . 

HANDLING  CONTAMINATION  .30000'  SPILLAGE  WEIGHT  LOSS  *00000 

HANDLING  CONTAMINATION  .50000'  SPILLAGE  WEIGHT  LOSS  *05000 

BLEND  1 PRODUCES  FOOD  ..ft  AT  TIME  17  MITH_S,P._ Q FROM  FOODS 
2 3 4 5 

HANDLING  CONTAMINATION  ,40000*  SPILLAGE  WEIGHT  LOSS  *02000 
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BLEND  2 PRODUCES  FOOD  8 AT  TIME  25  WITH  S.D,  0 FROM  FOODS 

6 ? - 

HANDU1I®1 CONTAMINATION  <36000*  uSPILLAGE’ WEIGHT  LOSS  .01000 

BLEND  3 PRODUCES  FOOD  U AT  TIME  31  WITH  S.D.  0 FROM  FOODS 

8 9 10 

HANDLING. CONTAMINATION  .30000*  SPILLAGE  WEIGHT  LOSS  .050Q0 

BLEND  4 PRODUCES  FOOD  13  AT  TIME  39  WITH  S.D.  0 FROM  FOODS 

11  12 

HANDLING  CONTAMINATION  .40000*  SPILLAGE  WEIGHT  LOSS  .O4OQ0 


FINAL  SCHEDULE 


FOOD 

1 

STARTS 

AT 

1=  1 

AND 

ENDS 

AT 

SPLIT 

1 

STOR, 

NO.= 

3 

FOOD 

2 

STARTS 

AT 

1=  1 

AND 

ENDS 

AT 

BLEND 

1 

STOR, 

NO.= 

2 

FOOD 

3 

starts 

AT 

1=  1 

AND 

ENDS 

AT 

BLEND 

1 

STOR. 

NO.= 

2 

FOOD 

4 

STARTS 

AT 

1=  13 

AND 

ENDS 

AT 

BLEND 

1 

STOR. 

NO  .= 

3 

FOOD 

5 

STARTS 

AT 

1=  1 

AND 

ENDS 

AT 

BLEND 

1 

STOR. 

NQ.= 

2 

FOOD 

6 

STARTS 

AT 

1=  17 

AND 

ENDS 

AT 

BLEND 

2 

STOR. 

NO.= 

4 

FOOD 

7 

STARTS 

AT 

1=  1 

AND 

ENDS 

AT 

BLEND 

2 

STOR. 

NO.= 

1 

FOOD 

e 

starts 

AT 

1=  25 

AND 

ENDS 

AT 

BLEND 

3 

STOR. 

NO.= 

4 

FOOD 

9 

STARTS 

AT 

1=  25 

AND 

ENDS 

AT 

blend 

3 

STOR. 

NO . = 

3 

FOOD 

10 

STARTS 

AT 

1=  13 

AND 

ENDS 

AT 

BLEND 

3 

STOR. 

NO . - 

2 

FOOD 

11 

STARTS 

AT 

1=  31 

AND 

ENDS 

AT 

BLEND 

4 

STOR. 

NO.= 

4 

FOOD 

12 

STARTS 

AT 

1=  35 

AND 

ENDS 

AT 

BLEND 

4: 

STOR. 

NO  . - 

4 

FOOD 

13 

STARTS 

AT 

1=  39 

AND 

IS  SERVED  AT  I 

97  STOR,  NO 

.=  3 

Js  1 

CHICKEN 

IC  = 

1 

IE=  13  1 

LS= 

1 LB= 

0 NATC^ 

1 

Ws  bO. 

0000  Yl: 

- 

5.000C 

I VI: 

: 6 

.0000  YA= 

22.3770 

K=  1 IACH=  4 YAC=100. 00000 

J=  2 SAUCE  (ROUX) 

IC=  1 IE=  17  LS=  0 LB=  1 NATC=  1 

W=  2.0000  YI=  5.0000  Vl=  6.0000  YA=  5*000Q 
K=  1 I ACH-  13  YAC=100. OOOOO 

J=  3 STOCK  AND  SEASONING 

IC=  1 IE=  17  LS=  0 LB-  1 NATC3  0 
W=  23.0000  YI=  5.0000  Vis  6.0000  YA~ 

J=  4 EXTRA  STOCK 

ICs  13  IE=  17  LS=  0 LB=  1 NATC=  0 
Ms  .0000  YI=  .0000  VI=  .0000  YA= 

J=  5 CHILLED  VEGETABLES 

1C=  1 IE=  17  LS=  0 LB=  1 NATC=  0 

Ws  8.0000  Ylr  5,0000  Vis  7.0000  YA= 

■ * : -i  . 

Js  6 MIXTURE  /w°.(A 

ICs  17  IE=  25  LS=  0 LB=  2 NATCs-  0 , , 

Ms  .0000  YI  = .0000  Vis  .0000  YAs 

r _ . 4 

J=  7 FROZEN  PEAS 

IC=  1 IE=  25  LSa  0 LBs  2 MK&l.  Q 

Ms  6.0000  YIS  “18.0000  Vis  5.0000  YA*  niB.OOOO 


5.0000 


21.9099 


5.0000 


95,8263 
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J=  e MIXTURE  BETA 

IC=  25  I£:  31  LS:  0 l.B=  3 NATCr  0 

*0000  YI:  .0000  VI:  .0000  YA:  114,0200 

J=  9 SLURRY 


1C=  25  IE:  31  LS:  0 LB:  3 NATC: 
W=  2.0000  YI:  10,0000  VI: 

K:  1 IACH:  27  YAC:100 . 00000 


1 

6.0000  YA: 


27.1119 


J=1°  OE-BONED*  COOKED  CHICKEN 

IC=  13  IE:  31  LS:  Q LB:  3 NATC:  1 

.0000  YI:  .0000  VI:  .0000  YA: 

K:  1 IACH:  27  YAC:  20.00000 


5.00Q0 


J=H  MIXTURE  GAMMA 

IC:  31  IE:  39  LS:  0 LB:  4 NATC:  0 

*l=  *0000  YI:  .0000  VI:  .0000  YA:  103.5712 

^=12  TOPPING 

IC:  35  IE:  39  LS:  0 LB:  4 NATC:  0 

5.0000  YI:  5.0000  VI:  6.0000  YA:  104.1283 

J=13  CHICKEN  POT  PIE 

IC:  39  IE:  97  LS:  0 LB:  0 NATC 

*=  .0000  YI:  .0000  VI: 

K:  1 IACH:  46  YAC:-18 . 00000 
K:  2 IACH:  88  YAC:l75 . 00000 


= 2 

,0000  YA:  19.2408 


1 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


FOOD  NUMBER  1 
Y:TEMP 
5*00000 
7.71663 
10.00856 
11.94218 
25*70872 
37.32306 
47*12167 
55*38841 
62*36277 
68.24680 
73*21094 
77.39901 
80.93234 
6.0000  RO: 


CHICKEN 
V:BAC . COUNT 
6.00000 
5.97024 


5.92052 
5.92616 
5.95663 
6.15091 
6.22938 
5.47188 
2.33034 
-4.96316 
-17.87334 
-37.35074 
-64.02174 
.0000  VE: 


A 

.00047 
.00112 
.00219 
.00371 
.05744 
.14503 
-.01001 
-1.46430 
-4.80428 
-9.83416 
-16.07734 
-23.02310 
-30123110 
-64.0217 


YA:AMBIENT  TEMP 
22.37698 
22.37698 
22.37698 
22.37698 
100.00000 
100.00000 
100.00000 
100.00000 
100.00000 
100.00000 
100.00000 
100.00000 
100.00000 


THIS  FOOD  ENDS  IN  SPLIT  1 AT  TIME*  1=  13.  INPUT  TO 

TEMP: Y:  80.932  BAC.  COUNT:V:  -64.022  WEI6HT:w: 

K AND  LINO(K)  10  21  31  40 

OUTPUT  FROM  SPLIT  IS  AS  FOLLOWS 

FOOD  NO.  4,  TIME:  13.  Y:  80.932.  V:  .216.  W: 

FOOD  NO.  10.  TIME:  13.  Y:  80,932.  V:  .289.  W: 


SPLIT  HAS 
60.000 


18.000 

35.918 


FOOD  NUMBER  2 
I Y:TEMP 

1 5.00000 

2 5.00000 

3 5.00000 


SAUCE  (ROUX) 


V:BAC.  COUNT  A 

6*00000  .00047 

6.00949  ,00047 

6.04303  .00047 


YA:AMBIENT  TEMP 
5.00000 
5.00000 
5.00000 


IAB 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


IAb 


1 

2 

3 
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4 5*00000 

6.01281 

.00047 

5.00000 

4 

5 5.00000 

6.13465 

>>00047 

5.00000 

5 

6 5.00000 

6.12334 

>00047 

5.00000 

6 

7 5.00000 

6.19452 

.00047 

5.00000 

7 

6 5*00000 

6.23755 

.00047 

5.00000 

8 

9 5.00000 

6.24146 

.00047 

5.00000 

9 

10  5.00000 

6.30309 

.00047 

5.00000 

10 

11  5.00000 

6.34361 

.00047 

5.00000 

11 

12  5*00000 

6.28498 

.00047 

5.00000 

12 

13  5.00000 

6.23829 

.00047 

5.00000 

13 

14  42.37959 

6.34732 

.14403 

100.00000 

14 

15  65.05145 

3.00684 

-6.81470 

100.00000 

15 

16  76*50264 

-13.11792 

— ? 25.73138 

1QQ. 00000 

16 

17  87.14315 

-49.13859 

-46.26512 

100.00000 

17 

RI  = 6.0000  RO= 

.0000  VE=  ■ 

-49.1366 

FOOD  NUMBER  3 

STOCK 

and  SEASONING 

I Y=TEMP 

V=BAC.  COUNT 

A 

YA=AMBIENT  temp 

IAB 

1 5.00000 

b. 09000 

.00047 

5.00000 

1 

2 5.00000 

5.99103 

.00047 

5.00000 

2 

3 5.00000 

6.02379 

.00047 

5.00000 

3 

4 5.00000 

6.05223 

.00047 

5.00000 

4 

5 5.00000 

6.06751 

.00047 

5.00000 

5 

b 5.00000 

6.14797 

.00047 

5.00000 

6 

7 5.0Q0QQ 

6.16938 

.00047 

5.00000 

7 

8 5*00000 

6.15180 

.00047 

5.00000 

8 

9 5.00000 

6.10543 

.00047 

5.00000 

9 

10  5.00000 

6.07439 

.00047 

5.00000 

10 

11  5.00000 

o. 04348 

.00047 

5.00000 

11 

12  5.00000 

6.12607 

.00047 

5.00000 

12 

13  5.00000 

6.18875 

.00047 

5.00000 

13 

14  5.00000 

6.17484 

.00047 

5.00000 

14 

15  5.00000 

6.24848 

.00047 

5.00000 

15 

16  5.00000 

6.23094 

.00047 

5.00000 

16 

17  5*00000 

6.16659 

.00047 

5.00000 

17 

RI=  6.0000  RO= 

6.0000  VE= 

6.1666 

MONREC  signal 

FOOD  NUMBER  4 

EXTRA 

STOCK 

I Y=TEMP 

V=BAC.  COUNT 

A 

ya=ambient  temp 

IAB 

1 80.93234 

*21616 

-.59223 

21.90990 

13 

2 57.70882 

-1.29140 

-2.31021 

21.90990 

14 

3 43*62304 

-2*33449 

.12893 

21.90990 

15 

4 35.07958 

-2.21477 

.13387 

21.90990 

16 

b 29.89771 

-2.03550 

.09285 

21.90990 

17 

RI=  .0000  RO= 

.4764  VE= 

-2.0355 

FOOD  NUMBER  5 

chilled  vegetables 

I YiTEMP 

VsBAC.  COUNT 

A 

YAzAMBIENT  temp 

IAB 

1 5.00000 

7.00000 

.00047 

5.00000 

1 

2 5.00000 

7.06141 

.00047 

5.00000 

2 

3 5.00000 

7.07714 

.00047 

5.00000 

3 

4 5.00000 

6.97111 

.00047 

5.00000 

4 

5 5.00000 

7.00327 

.00047 

5.00000 

5 

6 5.QQQQQ 

6*94106 

*00i*T 

5.80000 

6 

7 5.00000 

7.02143 

.00047 

5,00000 

7 

8 5*00000 

6.91424 

•00047 

5.00000 

8 
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9 

5.00000 

7.01987 

10 

5.00000 

7.08657 

11 

5*00000 

7.10434 

12 

5*00000 

7*09268 

13 

5.00000 

7.07536 

14 

5*00000 

7.13735 

15 

5*00000 

7*08740 

16 

5*00000 

7*12174 

17 

5*00000 

7.12076 

RI  = 

6.0000  RO: 

6.0000  VEr 

monrec 

SIGNAL 

,00047 

.00047 

.00047 

.00047 

.00047 

,00047 

.00047 

.00047 

.00047 

7.1208 


5.00000  9 
5.00000  10 
5.00000  11 
5.00000  12 
5.00000  13 
5.00000  14 
5.00000  15 
5.00000  16 
5.00000  17 


THIS  FOOD  ENDS  IN  BLEND  NO. 
FOOD  NO.  2,  TIME=  17,  Y = 
FOOD  NO.  3#  TIME:  17,  Y= 
FOOD  NO.  4,  TIME:  17,  Y= 
FOOD  NO.  5,  TIM£=  17,  Y= 
K AND  LINQ(K)  1 0 2 

OUTPUT  FROM  BLEND  HAS 

FOOD  NO.=  6,  TIME:  17,  Y= 


1.  INPUT  TO  BLEND  HAS 
87.143,  V:  -49.139#  WT.r 
5.000,  V:  6.167#  WT.: 

29.898,  V:  -2.035#  WT.r 

5.000#  V:  7.121#  WT.= 

3 1 4 1 

17.009#  V:  6.024#  WT.: 


2.000 

23.000 

18.000 

8.000 


49.344 


FOOD  NUMbER  6 
I Y=TEMP 

1 17.00873 

2 48.02103 

3 66.63095 

4 78. >3974 

5 85.15952 

6 89.35658 

7 91.90222 

8 93.44623 

9 94.38272 
RI=  3.1191  RO: 


FOOD  NUMBER  7 


I 

Y:TEMP 

1 

-18*00000 

2 

-18.00000 

3 

-18.00000 

4 

-18.0T000 

5 

-18*00000 

6 

-18.00000 

7 

-18.00000 

8 

-18*00000 

9 

-18.00000 

10 

-18.00000 

11 

-18.00000 

12 

-18.00000 

13 

-18.00000 

14 

-18*00000 

15 

-18.00000 

16 

-18*00000 

17 

-18*00000 

16 

-18.00000 

19 

-l&tQQOQQ 

20 

-18.00000 

21 

-18.00000 

MIXTURE 
V:BAC.  COUNT 
6.02421 
6.00319 
1.75707 
“14.74459 
-47.39148 
“94.23379 
-151.66308 
-216.30290 
-285.61029 
.0000  V£:  -285 

FROZEN  PI 
V:BAC.  COUNT 
5.00000 
4.97947 
4.84784 
4.89009 
4.98503 
5.00463 
4.91015 
4.95256 
4.93969 
4.95704 
5.00335 
5.06126 
5.01780 
5.03225 
5.02824 

5.11169 
5.13433 

5.11170 
5.0*629 
5.12079 
5*12691 


ALPHA 


A 

YA:AMBIENT  TEMP 

iab 

.01251 

95.82632 

17 

-.07861 

95.82632 

18 

-8.40681 

95.82632 

19 

“24.62128 

95.82632 

20 

“40.64919 

95.82632 

21 

“53.11461 

95.82632 

22 

“61.78862 

95.82632 

23 

“67.48243 

95.82632 

24 

-71.10028 

6103 

95.82632 

25 

AS 

A 

YA:AMBIENT  TEMP 

Iab 

.00000 

-18.00000 

1 

2 

3 

4 

5 

6 
7 
6 
9 

10 


.00000  -18.00000 

.00000  -18.00000 

.00000  -18.00000 

,00000  -18.00000 

.00000  -18.00000 

.♦00000  -18.00000 

.00000  -18.00000 

.00000  -18.00000 

,00000  -18.00000 

.00000  -18.00000  11 

.00000  -18.00000  12 

.00000  -18.00000  13 

.00000  -18.00000  14 

,00000  -18.00000  15 

.00000  -18.00000  16 

.00000  -18,00000  17 

.00000  -18.00000  18 

.00000  -18.00000  19 

.00000  -18.00000  20 

.00000  -18.00000  21 
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22 

-18*00000 

5.101Q9 

.00000 

-18,00000 

22 

23 

-18.00000 

5.14211 

.00000 

-18.00000 

23 

24 

-10.00000 

5*18777 

.00000 

-18.00000 

24 

25 

-18.00000 

5.15849 

.00000 

-18.00000 

25 

R I — 

6.0000  ROz 

6.0000  VE= 

5.1585 

monrec 

SIGNAL 

THIS  FOOD  ENOS  IN  BLEND  NO. 
FOOD  NO,  6,  TIME=  25,  Y = 
FOOD  NO.  7,  TIMEz  25,  Y= 
► and  LIND(K)  J i 2 
output  from  bleno  has 
^OOD  NO.=  S,  TIME;  25f  T= 


2.  input  to  bleno  has 

94.363,  Vz-285.610*  WT.z 
-18.000,  Vz  5.158,  WT.z 
0 3 0 4 1 

82.199,  Vs  3.433,  tfT.= 


49,344 

6.000 


55.092 


FOOD  NUMbER  8 
I YzTEMP 

1 82.19894 

2 92.68970 

3 99.72187 

4 104 , 43567 

5 107.59543 

6 109.71348 

7 111.13325 

RI  = 3.0000  ROz 


MIXTURE  BETA 
VzBAC.  COUNT 
3.43263 


-29.18136 
-106.65348 
-215.01946 
-342.60241 
-485.96942 
-640.52156 
,0000 


A 

-.63687 

-64.65088 

-94.21347 

-118.36759 

-136.67193 

-149.94077 

-159.29894 


VE=  -640.5216 


YAzAMBIENT  TEMP  Iab 

114.01999  25 

114.01999  26 

114.01999  27 

114.01999  28 

114.01999  29 

114.01999  30 

114.01999  31 


FOOD  NUMBER  9 SLURRY 


I 

YzTEMP 

VzBAC.  COUNT 

A 

YAzAMBIENT  TEMP 

iab 

1 

10.00000 

6.00000 

.00218 

27.11190 

25 

2 

15 » b4l45 

6.00523 

.00923 

27.11190 

26 

3 

19.42303 

6.00461 

.02050 

27.11190 

27 

4 

45  * 98764 

6.12632 

,05490 

100,00000 

28 

5 

63.79443 

3.26576 

-5.81801 

100.00000 

29 

6 

75.73068 

-9.72106 

-20.05820 

100.00000 

30 

7 

83.73179 

-38.23226 

-36.90062 

100.00000 

31 

RI  = 

6.0000  ROz 

,0000  VEz  -38.2323 

FOOD  NUMBER  10 

de-boned, 

COOKED  CHICKEN 

I 

YzTEMP 

VzBAC.  COUNT 

A 

YAzAMBIENT  TEMP 

iab 

1 

60.93234 

.28885 

-.59223 

5.00000 

13 

2 

69.06144 

-5.33824 

-10.72315 

5.00000 

14 

3 

59*04639 

-12.16957 

-2.91085 

5.00000 

15 

4 

50.59703 

-13.85776 

-.37718 

5.00000 

16 

5 

43.46661 

-13.95789 

.13151 

5.00000 

17 

6 

37.45462 

-13.89031 

,14550 

5.00000 

18 

7 

32.38082 

-13.81857 

.11418 

5.00000 

19 

8 

28.10023 

-13.76476 

.07710 

5.0000C 

20 

9 

24.48885 

-13.69590 

.04842 

5.00000 

21 

10 

21.44206 

-13.60814 

.02972 

5.00000 

22 

11 

18  .£7159 

-13.63677 

.01840 

5.00000 

23 

12 

16.70297 

-13.59368 

.01170 

5.00000 

24 

13 

14.87338 

-13.51709 

.00772 

5.00000 

25 

14 

13.32983 

-13.50370 

.00530 

5.00000 

26 

15 

12.02758 

-13*50436 

.00379 

5.00000 

27 

16 

13.27395 

-13.49092 

.00523 

20.00000 

20 

J.7 

TltiW- 

-IX*  41179 

20.00000 

29 

16 

15*21260 

-13.44995 

* 00036 

20.00000 

30 

19 

15.96104 

-13.49983 

»Q0§92 

20.00000 

31 

41 


RI= 


.0000  ROs  1,8000  V£s  -13.4998 


THIS  FOOD  ENOS  IN  BLEND  NO.  3.  INPUT  TO  BLEND  HAS 


FOOD 

NO.  8.  TIME= 

31,  Ys  111.133, 

Vs-640.522, 

WT.s  55.092 

FOOD 

NO.  9.  TIME= 

31?  Y=  83.732 i 

V“  -38.232' 

wt.s  2.000 

FOOD 

NO.  10*  TIMES 

31,  Ys  15.961* 

Vs  -13.500, 

WT.s  35.918 

K ANO  LINQ(K)  1 0 

2 1 3 

.1  4 

1.  .. 

OUTPUT 

FROM  BLEND  HAS 

FOOD 

NO.sil,  TIMES 

31*  Y=  73.791* 

Vs  .191, 

WT.s  89,711 

FOOD  NUMBER  11 

MIXTURE  < 

GAMMA 

I 

YsTEMP 

VsBAC.  COUNT 

A YAsAMBIENT  TEMP 

iab 

1 

73*79101 

,19133  ..  . 

-.35583 

103.57116 

31 

2 

81.50949 

-15.74546 

-31.53442 

103.57116 

32 

3 

87.22747 

-54.76420 

-46.51463 

103,57116 

33 

4 

91.46346 

-108.11862 

-60.23105 

103.57116 

34 

5 

94.60156 

-174.10405 

-71.96387 

103,57116 

35 

6 

96*92632 

-250.90991 

-81.57274 

103.57116 

36 

7 

98.64854 

-336.32293 

-89.21661, 

103.57116 

37 

6 

99.92439 

-428.51825 

-95.17672 

103.57116 

38 

9 

100.86957 

-525.97297 

-99.75896 

103.57116  . 

39 

RI  = 

.6000  ROs 

.0000  VEs  -525 

.9730 

FOOD  NUMBER  12 

topping 

I 

YsTEMP 

VsBAC.  COUNT 

A YAsAMBIENT  TEMP 

IAB 

1 

5*00000 

6.00000 

.00047 

104.12827 

35 

2 

20.49724 

6.04972 

.02510 

104.12827 

36 

3 

33.57171 

6.11109 

.12352 

104.12827 

37 

4 

44.60218 

6.22913 

.10682 

104.12827 

38 

5 

53.90820 

5.72607 

-1.04184 

104.12827 

39 

RI  = 

6.QQ0Q  ROs 

.0000  VEs  5 

.7261 

THIS  FOOD  ENDS  IN  BLEND  NO.  4.  INPUT  TO  BLEND  HAS 

FOOD 

NO,  11*  TIMES 

39*  Ys  100.870. 

V=-525.973, 

WT.s  89.711 

FOOD 

NO.  12.  TIMES 

39*  Ys  53,908, 

Vs  5.726* 

WT.s  5.000 

K AND  1 

LIND(K)  1 0 

2 0 3 

1 4 

1 

OUTPUT 

FROM  BLEND  HAS 

food 

NO. =13.  TIMES 

39*  Ys  98.390, 

Vs  3.307, 

WT.s  86.701 

FOOD  NUMBER  13 

chicken 

POT  PIE 

I 

YsTEMP 

VsBAC.  COUNT 

A YAsAMBIENT  TEMP 

IAB 

1 

98.39039 

3.30671 

-1.36038 

19.24077 

39 

2 . 

84.04300 

-16.23512 

-37.69725 

19.24077 

40 

3 

72.29635 

-42.55360 

-14.77327 

19,24077 

41 

4 

62.67900 

-52.41697 

-5.01744 

19.24077 

42 

5 

54.80499 

-55.55599 

-1.28757 

19,24077 

43 

6 

48.35829 

-56.26345 

-.10868 

19.24077 

44 

7 

43*08018 

-56.23576 

.13705 

19.24077 

45 

„ a 

. 38.75883 

-56,01516 

•14892 

19.24077 

46 

9 

28.47020 

-55.85092 

.08031 

-18.00000 

47 

10 

20*04656 

-55.79982 

, .02309 

-18.00000 

48 

n 

13.14990 

-55.81018 

.00506 

-18.00000 

49 

12 

7.50338 

-55.82317 

.00105 

-18.00000 

50 

13 

2.88041 

-55.77293 

,00023 

-18.00000 

51 

- 14- 

..  -.90457- 

-5 5.7952 Q 

- . *oihlqq 

.^la.Qoaoo 

52 

15 

-4.00345 

-55.77630 

.00000 

-18,00000 

53 

16 

-6*54059 

-55,75830 

.00000 

-18.00000 

54 

42 


17 

-8.61783 

-55.75044 

.00000 

-18,00000 

55 

18 

-10.31853 

-55.68500 

,00000 

-18.00000 

56 

19 

-11.71094 

-55.63998 

,00000 

-18.00000 

57 

20 

-12.85096 

-55.69514 

.00000 

-16.00000 

58 

21 

-13.78432 

-55.68975 

.00000 

-18.00000 

59 

22 

-14.54849 

-55.59563 

.00000 

-18.00000 

60 

23 

-15.17414 

-55.54278  _ 

.00000 

-18.00000 

61 

24 

-15.68639 

-55.61711 

.00000 

-18.00000 

62 

2b 

-16.10577 

-55.77279 

,00000 

-18.00000 

63 

26 

-16*44914 

-55.74161 

.00000 

-.8.00000 

64 

27 

-16.73026 

-55.79869 

.00000 

-16.00000 

65 

26 

-16.96043 

-55.69548 

.00000 

-18.00000 

66 

29 

-17.14887 

-55*61345 

•QQQQQ 

rlB. 00000 

67 

30 

-17.30315 

-55.58629 

.00000 

-18.00000 

68 

31 

-17.42947 

-55.62767 

.00000 

-18.00000 

69 

32 

-17.53289 

-55.70799 

.00000 

-18.00000 

70 

33 

-17.61756 

-55*66743 

.00000 

-18.00000 

71 

34 

-17.68689 

-55.6 c.  729 

.00000 

-18.00000 

72 

35 

-17.74364 

-55.69900 

.00000 

-18.00000 

73 

3b 

-17.79011 

-55.78619 

.00000 

-18.00000 

74 

37 

-17.82816 

-55.67530 

.00000 

-18.00000 

75 

38 

-17.65931 

-55.62170 

.00000 

-18.00000 

76 

39 

-17.88481 

-55.58746 

.00000 

-18.00000 

77 

40 

-17.90569 

-55.58092 

.00000 

-18.00000 

78 

41 

-17.92279 

-55*48587 

.00000 

-18.00000 

79 

42 

-17.93676 

-55.52533 

.00000 

-18.00000 

80 

43 

-17.94824 

-55.57060 

.00000 

-18.00000 

81 

44 

-17.95762 

-55.48525 

.00000 

-18.00000 

82 

45 

-17.96531 

-55.47280 

.00000 

-18.0000C 

83 

46 

-17.97160 

-55.45072 

.00000 

-18.00000 

84 

47 

-17.97674 

-55.48998 

.00000 

-18.00000 

85 

46 

-17.98096 

-55 . 56344 

.00000 

-18.00000 

86 

49 

-17.98441 

-55.47215 

,00000 

-18.00000 

87 

50 

-17.98724 

-55.43489 

.00000 

-18.00000 

88 

51 

16.99541 

-55.56133 

.01247 

175.00000 

89 

52 

45.63679 

-55.54091 

.07066 

175.00000 

90 

53 

69.08636 

-60.86979 

-10.75114 

175.QQQ00 

91 

54 

88.26525 

-91.17761 

-49.72088 

175.00000 

92 

55 

104.00396 

-173.92798 

-116.00137 

175.00000 

93 

56 

116.87336 

-332.54395 

-201.09789 

175.00000 

94 

57 

127.40993 

-580.95062 

-295.70990 

175.00000 

95 

58 

136.03655 

-924.92157 

-392.23589 

175.00000 

96 

59 

143.09942 

-1363.7011Q 

-485.33123 

175.00000 

97 

THIS 

FOOD  IS  SERVED 

AT  TIME  97 

R I- 

.0000  R0= 

.0000  VE--1363 

,.7011 

MAX. 

COUNT  IS  7,1373  IN  FOOD  NO. 

5 AT  TIME  14 

NSIG=1  NSIGD=0  NTRUrO 
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FINAL  schedule 


FOOD 

1 

STARTS 

AT 

Ir 

1 

AND 

FOOD 

2 

STARTS 

AT 

Ir 

1 

AND 

FOOD 

3 

STARTS 

AT 

Ir 

1 

AND 

FOOD 

4 

STARTS 

AT 

Ir 

13 

ANO 

FOOD 

5 

starts 

AT 

Ir 

1 

AND 

FOOD 

6 

STARTS 

AT 

Ir 

17 

AND 

FOOD 

7 

STARTS 

AT 

Ir 

1 

ANO 

FOOD 

6 

STARTS 

AT 

Ir 

25 

AND 

FOOD 

9 

STARTS 

AT 

Ir 

25 

AND 

food 

10 

STARTS 

at 

Ir 

13 

AND 

FOOD 

11 

starts 

AT 

Ir 

31 

AND 

FOOD 

12 

STARTS 

AT 

Ir 

35 

AND 

FOOD 

13 

STARTS 

AT 

Ir 

39 

AND 

ENDS  AT  SPLIT  1 STOR.  NO.r  3 
ENDS  AT  BLEND  1 STOR.  NO,:  2 

ENDS  AT  BLEND  1 STOR.  NO.r  2 

ENDS  AT  BLEND  l.STOR.  NO.r  3 
ENDS  AT  BLEND  1 STOR.  NO.r  2 

ENDS  AT  BLEND  2 STOR.  NO.r  4 

ENDS  AT  BLEND  2 STOR.  NO.r  1 

ENDS  AT  BLEND  3 STOR.  NO.r  4 

ENDS  AT  BLEND  3 STOR.  NO.r  3 

ENDS  AT  BLEND  3 STOR.  N0.r  2 

ENDS  AT  BLEND  4 STOR.  NO.r  4 

ENDS  AT  BLEND  4 STOR.  NO.r  4 

IS  SERVED  AT  Zr  97  STOR.  NO.=  3 


J=  1 CHICKEN 

lCr  1 IEr  13  LSr  1 LBr  0 NATC=  1 

Wr  60.0000  Ylr  5,0000  Vlr  6.0000  YAr 

Kr  i IACHr  4 YACrlOO . 00000 


21.6280 


Jr  2 
IC  = 
Mr 

Kr  1 


SAUCE  (ROUX) 

1 IEr  17  LSr  0 LBr  1 NATCr  1 
2.0000  Ylr  5.0000  Vlr  6.0000  YAr 
IACHr  13  YACrlOO. 00000 


5.0000 


J=  3 STOCK  AND  SEASONING 

IC=  1 IEr  17  LSr  0 LBr  1 NATCr  0 
Mr  23.0000  Ylr  5.0000  Vlr  6.0000  YAr  5.0000 

J=  4 EXTRA  STOCK 

ICr  13  IEr  17  LSr  0 LBr  1 NATCr  0 

Mr  18.0000  Ylr  80.2978  Vlr  .HU  YAr  21.9506 

J=  5 CHILLED  VEGETABLES 

ICr  1 IEr  17  LSr  0 LBr  1 NATCr  Q 

Mr  8.0000  Ylr  5.0000  Vlr  7.0000  YAr  5.0000 

J=  6 mixture  alpha 

ICr  17  IEr  25  LSr  0 LBr  2 NATCr  0 

*=  ‘♦9.8154  Ylr  13.7372  Vlr  5.9348  YAr  98.7461 

J=  7 . FROZEN  PEAS 

IC=  1 IEr  25  LSr  0 LBr  2 NATCr  0 . 

Mr  6.0000  Ylr  -18.0000  Vlr  5.0000  YAr  -18.0000 

J=  8 MIXTURE  BETA 

IC=  25  IEr  31  LSr  0 LBr  3 NATCr  0 

Mr  55.1513  Ylr  100.5085  Vlr  2.9985  YAr  102.2474 


J=  9 SLURRY 

K=  25  IEr  31  LSr  Q LBr  3 NATCr  1 
Mr  2.0000  Ylr  10,0000  Vlr  6.0000  YAr 
Kr  1 IACHr  27  YACrlOO. QQQOO 


21.6481 


J=10  DE-BONED,  COOKED  CHICKEN 
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IC=  13  IE=  31  L§=  0 LB=  3 NATC:  1 

m=  33.4602  YI:  11.9689  Vl=  .3913  YA=  5.0000 

K=  1 IACH=  27  YAC=  20  00000 

J=ll  MIXTURE  GAMMA 

IC  = 31  IE=  39  LS:  0 LB:  4 NATC:  0 

W=  88.5702  Yl  = 71.9381  VI=  .0524  YA:  100,6015 

0=12  TOPPING 

1C=  3b  IE=  39  LS=  0 LB:  4 NATO:  0 

W = 5.0000  YI:  5,0000  VI:  6.0000  YA:  97.8588 


0=13 

CHICKEN  POT  PIE 

IC=  39 

IE:  97  LS:  0 

LB:  0 NATC:  2 

ft:  89 

.8401  YI:  -17 

.9916  VI=  3.4269 

YA:  23.8291 

K=  1 

IACH:  45  YAC=- 

18.00000 

K:  2 

IACH:  93  YAC:175 . 00000 

FOOD  NUMBER  1 

chicken 

I 

Y=TEMP 

V: BAC.  COUNT 

A YA:AMBIENT  TEMP 

IAB 

1 

5*00000 

6.00000 

.00047  21.62800 

1 

2 

7.59954 

6.00614 

.00108  21.62800 

2 

3 

9*79268 

5.93523 

.00206  21.62800 

3 

4 

11.64296 

5.97603 

.00343  21.62800 

4 

5 

25*45627 

6.01859 

.05551  100.00000 

5 

6 

37*11008 

6.15310 

.14423  100.00000 

6 

7 

46*94199 

6.29487 

.00180  100,00000 

7 

8 

55*23682 

5.54775 

-1,41706  100.00000 

8 

9 

62*23488 

2.48154 

-4.71977  100.00000 

9 

10 

68.13890 

-4.70745 

-9.72018  100.00000 

10 

11 

73*11991 

-17.55583 

-15.94427  100.00000 

11 

12 

77*32221 

-36.99811 

-22.88064  100.00000 

12 

13 

80*86755 

-63.42109 

-30.08706  100,00000 

13 

RI  = 

6.0000  RO: 

.0000  VE=  - 

63. 

4211 

THIS 

FOOD  ENDS  IN  SPLIT  1 AT  TIME 

, I 

= 13.  INPUT  TO  SPLIT  HAS 

temp 

= Y=  80.668 

BAC.  COUNT=V= 

-63 

.421  WEIGHT:*:  60,000 

K ANO  1 

LIND(K)  1 0 

2 1 

3 

1 4 0 

OUTPUT 

FROM  SPLIT  IS  AS  FOLLOWS 

food 

NO.  4,  TIME 

: 13*  Y=  80.868, 

V:  .095*  W=  18.000 

FOOD 

NO,  ig»  TIME 

= 13*  Y=  80,868 r 

V:  ,000*  W=  34.739 

FOOD  NUMBER  2 

SAUCE 

(ROUX) 

I 

Y=TEMP 

V:BAC.  COUNT 

A YA:AMBIENT  TEMP 

I AB 

1 

5.00000 

6*00000 

,00047  5,00000 

1 

2 

5*00000 

5.99298 

.00047  5.00000 

2 

3 

5*00000 

5.93149 

.00047  5.00000 

3 

4 

5.00000 

5.95563 

.00047  5,00000 

4 

5 

5t0Q000 

5,97959 

.00047  5.00000 

5 

6 

5.00000 

6.03102 

.00047  5.00000 

6 

7 

5*00000 

6.05524 

,00047  5.00000 

7 

6 

5*00000 

6.04468 

.00047  5,00000 

e 

9 

5*00000 

6.00616 

,00047  5,00000 

9 

10 

5.09000 
. youfl. 

6.02130 

,00047  5.00000 

10 

11 

6,01941 

- 5*000011 

,00047  5.00000 

U 

12 

5.00000 

5.90511 

12 

13 

5*00000 

5.50074 

,00047  5,00000 

13 

45 


14 

42.37959 

5.90162 

15 

65*05145 

2.53755 

16 

78*80264 

-13*68201 

17 

87*14315 

-49.75846 

RI  = 

6.0000  R0= 

.0000  VE= 

FOOD  NUMBER  3 

STOCK 

I 

Y=T£MP 

V=8AC.  COUNT 

1 

5.00000 

6.00000 

2 

5.00000 

5.99296 

3 

5*00000 

5.98290 

4 

5.00000 

5.97366 

5 

5*00000 

5.94438 

6 

5*00000 

5.88292 

7 

5*00000 

5.93371 

8 

5*00000 

5.88202 

9 

5*00000 

5.89114 

10 

5*00000 

5.91858 

11 

5.00000 

5.93891 

12 

5.00000 

5.92722 

13 

5 * OOOoO 

5.94820 

14 

5.00000 

6.00450 

15 

5*00000 

6.03992 

16 

5*00000 

5.97079 

17 

5*00000 

5.98977 

RI  = 

6.0000  RO= 

6.0000  VE= 

monrec 

SIGNAL 

.14403 

100.00000 

14 

-6.81470 

100.00000 

15 

-25.73138 

100.00000 

16 

-46.26512 

100.00000 

17 

-49.7585 

AND  SEASONING 

A YArAMBlENT  TEMP  IAB 

.00047  5.00000  1 

.00047  5.00000  2 

.00047  5.00000  3 

.00047  5.00000  4 

.00047  5.00000  5 

.00047  5.00000  6 

.00047  5.00000  7 

.00047  5.00000  8 

.00047  5.00000  9 

.00047  5.00000  10 

.00047  5.00000  U 

.00047  5.00000  12 

.00047  5.00000  13 

.00047  5.00000  14 

.00047  5.00000  15 

.00047  5.00000  16 

,00047  5.00000  17 

5.9898 


FOOD  NUMBER  4 
I YsTEMP 

1 80.86755 

2 57.68553 

3 43.62492 

4 35.09673 

5 29.92412 
RI=  .0000  R0= 


EXTRA 
V=BAC.  COUNT 

STOCK 

A 

YArAMBlENT  TEMP 

IAB 

•09506 

-.58988 

21.95057 

13 

-1.25051 

-2.30052 

21.95057 

14 

-2.34451 

.12890 

21,95057 

15 

-2.24988 

.13397 

21.95057 

16 

-2.16457 

.09308 

21.95057 

17 

.4780  VE= 

-2.1646 

FOOD  NUMBER  5 CHILLED  VEGETABLES 


I 

Y=TEMP 

V=BAC.  COUNT 

1 

5.00000 

7.00000 

2 

5*00000 

7.05232 

3 

5*00000 

7.06688 

4 

5*00000 

7.09419 

5 

5*00000 

7.10400 

6 

5*00000 

7.08599 

7 

5*00000 

7.03965 

8 

5*00000 

7.14487 

9 

5*00000 

7.25272 

10 

5.Q0000 

7.32622 

11 

5.00000 

7.37216 

12 

5*00000 

7.40108 

13 

5*00000 

7.34428 

14 

5*00000 

7.33122 

15 

5*00000 

7.21676 

16 

. 5*128000 

-7j 17244 

17 

5*00000 

7.09677 

= 

6.0000  RO= 

6.0000  VE= 

A 

YArAMBlENT  TEMP 

IAB 

.00047 

5.00000 

1 

.00047 

5.00000 

2 

.00047 

5.00000 

3 

. .00047 

5.00000 

4 

.00047 

5.00000 

5 

,00047 

5.00000 

6 

.00047 

5.00000 

7 

.00047 

5.00000 

8 

.00047 

5.00000 

9 

•Q0047 

5.00000 

10 

.00047 

5,00000 

11 

. .00047 

5.00000 

12 

.00047 

5.00000 

13 

.00047 

5.00000 

14 

.00047 

5.00000 

15 

_ _•  Q0047 

4*8110(10 

16 

•00047 

7.0968 

5.00000 

17 
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monrec  signal 


THIS  FOOD  ENDS  IN  BLEND  NO.  i.  INPUT  TO  BLEND  HAS 


FOOD 

NO. 

2* 

TIME: 

17  * 

Y: 

87.143* 

V: 

-49.758* 

WT.r 

2.000 

FOOD 

NO. 

3 * 

TIME: 

17 * 

Y= 

S.OQO* 

V: 

5.990* 

WT.: 

23,000 

FOOD 

NO. 

4 * 

TIME: 

17  * 

Y= 

29.924* 

V: 

-2.165* 

WT.r 

18.000 

FOOD 

NO. 

5* 

TIME: 

17* 

Y= 

5.QHQ* 

V: 

7.097* 

WT.: 

8.000 

K ANO  LIND(K) 

1 0 

2 

1 3 

1 

4 

1 

OUTPUT 

FROM 

BLENO  HAS 

FOOD 

NO.= 

6 * 

TIME: 

17* 

Y = 

17.018, 

V: 

6.075* 

WT.r 

50.659 

BAD  SMELL  DETECTED 


FOOD  NUMBER  6 

MIXTURE 

ALPHA 

I 

Y:TEMP 

V:BAC.  COUNT 

A 

YA=AMBIENT  TEMP 

IAB 

1 

17.01805 

6.07520 

,01253 

98.74610 

17 

2 

49.17553 

6.04787 

-.19221 

98.74610 

18 

3 

68*68003 

.04899 

-10,30063 

98.74610 

19 

4 

80.51011 

-18.96976 

-29.30058 

98.74610 

20 

5 

87.68541 

-57.55078 

-47.88528 

90.74610 

21 

6 

92.03745 

-112.62434 

-62.27403 

98.74610 

22 

7 

94.67710 

-179.84624 

-72.26360 

98.74610 

23 

8 

96.27813 

-255.39060 

-78.81282 

90.74610 

24 

9 

97.24920 

-336.20288 

-82.97123 

98.74610 

25 

RI  = 

3.1195  RO: 

.0000  VE:  -336 

.2029 

FOOD  NUMBER  7 FROZEN  PEAS 


I 

Y:TEMP 

V:BAC • COUNT 

A 

YA:AMBIENT  TEMP 

IAI 

1 

-18.00000 

5.00000 

.00000 

-18.00000 

1 

2 

-18.00000 

5.04637 

.00000 

-18.00000 

2 

3 

-18.00000 

5.07207 

.00000 

-18.00000 

3 

4 

-18.00000 

5.07231 

.00000 

-18.00000 

4 

5 

-18.00000 

5.10884 

.00000 

-18.00000 

5 

6 

-18.00000 

5.19049 

,00000 

-18.00000 

6 

7 

-18.00000 

5.24387 

.00000 

-18.00000 

7 

8 

-18.00000 

5.20385 

,00000 

-18.00000 

8 

9 

-18.00000 

5.21094 

.00000 

-18.00000 

9 

10 

-18.00000 

5.19269 

,00000 

-18.00000 

10 

11 

-18.00000 

5.13368 

,00000 

-18.00000 

11 

12 

-18.00000 

5*06178 

.00000 

-18.00000 

12 

13 

-18.00000 

5.00474 

.00000 

-18.00000 

13 

14 

-18.00000 

4.93552 

,00000 

-18.00000 

14 

15 

-le.oopoo 

4.96833 

.00000 

-18,00000 

15 

16 

-10.QQQQQ 

5.02056 

. . QOOOQ 

-la.ocooo 

16 

17 

-18.00000 

5.00300 

.00000 

-18,00000 

17 

18 

-18.00000 

4.91601 

.00000 

-18.00000 

Id 

19 

-18*00000 

4.87908 

.00000 

-18.00000 

19 

20 

-18.00000 

4.94492 

.00000 

-18.00000 

20 

21 

-18.00000 

5.03349 

.00000 

-18.00000 

21 

22 

-18.QQQQQ 

5.02569 

.QQOaO 

-18.00000 

22 

23 

-18.00000 

5.08601 

.00000 

-18.00000 

23 

24 

-18.00000 

5.07235 

.00000 

-18.00000 

24 

25 

-18.00000 

5.07710 

.00000 

-18.00000 

2b 

RI  = 

6.0000  RO: 

6.0000  VE:  5. 

0771 

MONREC 

SIGNAL 

THIS  FOOD  ENOS  IN  BLEND  NO.  2.  INPUT 

TO  BLEND 

HAS 

FOOD 

NO,  6*  TIME: 

25*  Y=  97.249 » 

V:-336.203 

* WT.:  50.659 

47 


FOOD  NO,  7#  TIME:  25#  Y=  -18.000#  V: 
K AND  LIND(K)  11  20  30 
OUTPUT  FROM  BUENO  HAS 

food  no.:  a#  time:  25#  y=  85.045#  v= 


5,077»  *T,:  6,000 

4 1 

3.003#  WT,=  56.370 


1 
1 

2 

3 

4 

5 

6 
7 


RI  = 


FOOD  NUMBER  8 MIXTURE  BETA 


Y=TEMP 

V=BAC.  COUNT 

A 

YA=AMBIENT  TEMP 

IAB 

85 « 04467 

3.00321 

-.74886 

102.24743 

25 

90.71607 

-26.17553 

-57.63839 

102.24743 

26 

94.51773 

-90.82529 

-71.63224 

102.24743 

27 

97.06606 

-167.74942 

-82.17606 

102.24743 

28 

98.77425 

-253.75765 

-89.79252 

102,24743 

29 

99.91929 

-346.26971 

-95.15236 

102.24743 

30 

100.68683 

-443.23591 

-98.86183 

102.24743 

31 

0000  RO: 

.0000  VE=  -443 

.2359 

1 
1 

2 

3 

4 


5 

6 
7 


R I - 


1 NUMBER  9 
YrTEMP 

SLURRY 

VrBAC.  COUNT 

A 

YAzAMBIENT  temp 

IAB 

10.UO0QQ 

6.00000 

.00218 

21.64810 

25 

13.84014 

6.04010 

.00601 

21.64810 

26 

16.41427 

6.00773 

.01098 

21.64810 

27 

43.97081 

6.00166 

.12226 

100.00000 

28 

62.44251 

3.57010 

-4.85747 

100.00000 

29 

74.82446 

-8.11143 

-18.55928 

100.00000 

30 

83.12433 

-35.05844 

-35.37801 

100.00000 

31 

0000  RO: 

.0000  VE:  -35. 

0564 

FOOD  NUMBER  10 


I 

YrTEMP 

1 

80.66755 

2 

69.00678 

3 

59.00027 

4 

50.55813 

5 

43.43579 

6 

37.42692 

7 

32.35745 

8 

28*08052 

9 

24.47 222 

10 

21.42803 

11 

18.85975 

12 

16.69298 

13 

14.86496 

14 

13.32272 

15 

12.02159 

16 

13.26889 

17 

14.32120 

16 

15.20900 

19 

15.95800 

: 

.0000  RO: 

DE-BONED#  COOKED  CHICKEN 


:BAC . COUNT 

A 

YArAMBlENT  TEMP 

.00047 

-.58988 

5.00000 

-5.63722 

-10.66191 

5.00000 

-12.34975 

-2.88867 

5.00000 

-13.97480 

-.37139 

5.00000 

-14.10730 

.13203 

5.00000 

-13.91529 

.14540 

5.00000 

-13.71923 

.11399 

5.00000 

-13.62988 

.07693 

5.00000 

-13.51782 

.04830 

5.00000 

-13.50703 

.02965 

5.00000 

-13.38341 

.01836 

5.00000 

-13.39185 

.01168 

5.00000 

-13.42661 

.00770 

5.00000 

-13.43201 

.00529 

5.00000 

-13.46702 

.00379 

5.00000 

-13.49615 

.00522 

20.00000 

-13.50061 

.00676 

20.00000 

-13.51165 

.00835 

20.00000 

-13.52981 

.00991 

20.00000 

1,8000  VE=  -13.5298 


THIS  FOOD  ENDS  IN  BLEND  NO. 
FOOD  NO.  8 # TIME=  3l»  Y = 
FOOD  NO.  9»  TIME=  31#  Y= 
FOOD  NO.  10.  TIME:  31 » Y= 
K AND  LIND(K)  10  2 

OUTPUT  FROM  BLENO  HAS 

FOOD  NQ.=11»  TIME:  31#  Y= 


3.  INPUT  TO  BLEND  HAS 


100.687# 

V:-443.236»  WT .: 

56.370 

83 • 124  § 

V:  -35.058#  WT.= 

2.000 

15.958# 

V:  -13.530#  WT.r 

34.739 

1 3 

1 4 1 

68.698# 

V:  .096#  WT.r 

66.621 

IAB 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
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FOOD  NUMBER  11 

MIXTURE 

GAMMA 

I 

Y=TEMP 

VtBAC.  COUNT 

A YA 

sAMBIENT  temp 

IAB 

1 

68.69754 

.09593 

-.21483 

100.60149 

31 

2 

76.96646 

-11.21416 

-22.22831 

100,60149 

32 

3 

83.09223 

-39.99038  .. 

-35.29872 

100.60149. 

33 

4 

87.63031 

-81.53288 

-47.71895 

100,60149 

34 

5 

90.99220. 

r 134. 64308 

jr.5fi*58744 

100,60149 

.35 

6 

93.48275 

-197.79175 

-67.62117 

100.60149 

36 

7 

95.32780 

-269.01596 

-74,87982 

100.60149 

37 

8 

96.69464 

-346.69839 

-80,57905 

100.60149 

38 

9 

97.70723 

-429.47360 

-84.98229 

100.60149 

39 

RI  = 

.6000  RO= 

,0000  V E=  -429 

.4736 

FOOD  NUMBER  12 

topping” 

in  . , ~ - 

- 

X 

Y=TEMP 

V= BAC.  COUNT 

A YA 

sAMBIENT  temp 

IAB 

1 

5.00000 

6.00000 

.00047 

97.85880 

35 

2 

19.51710 

6.10702 

.02008 

97.85880 

36 

3 

31.76466 

6.20258 

,10906 

97.85880 

37 

4 

42j»  09750 

6.28747 

_•  14591  .. 

97,8588.0 

38 

5 

50.81495 

6.21084 

-.41037 

97.85880 

39 

Rl  = 

6.0000  RO= 

, 0Q00  VE=  6 

.2108 

THIS  FOOD  ENDS  IN  BLEND  NO.  4.  INPUT  TO  BLEND  HAS 

FOOD 

NO.  U*  TIME= 

39.  Y=  97.707* 

V=-429.474* 

WT.=  88.621 

FOOD 

NO,  12  i.  TIME? 

J9f  1=  5Q.,815j.  Y.5 8*211*  . 

WT.=  5.000  .... 

K AND 

LINO (K)  1 0 

2 0 3 

1 4 1 

OUTPUT 

FROM  BLEND  HAS 

FOOD 

NO. =13.  TIME= 

39*  Y=  95.203* 

V=  4.058* 

WT ,=  87.776 

FOOD  NUMbER  13 

chicken 

POT  PIE 

I 

Y=TEMP 

VrBAC.  COUNT 

w.  A ....  YA 

=AMBIENT  TEMP 

IAB 

1 

95.20285 

4.05788 

-1.19998 

23.82914 

39 

2 

82.26499 

-13.23796 

-33.29614 

23.82914 

40 

3 

71.67237 

-36.84569 

-13.92490 

23.82914 

41 

4 

62.99986 

-46.39201 

-5.23988 

23,82914 

42 

5 

55.89942 

-49.76034 

-1.63052 

23.82914 

43 

6. 

5Q.PJ36Q6_ 

-50.78745 

.....  -.30454 

23,82914 

44 

7 

45.32649 

-50.66390 

.08305 

23,82914 

45 

8 ' 

33.84734 

-50.60722 

.12555 

-18,00000 

46 

9 

24.44902 

-50.54232 

.04814 

-18,00000 

47 

10 

16.75431 

-50.48440 

.. .01163 

-18.00000 

48 

11 

10.45443 

-50.46622 

.00248 

-ia. ooooo 

49 

12 

5.29651 

-50,38521 

00 052. 

-18,00000 

50 

13 

1.07357 

-50.33064 

.00012 

-18.00000 

51 

14 

-2.38388 

-50.32083  _ 

.00000 

-la.OOQOQ 

52 

15 

-5.21460 

-50.35106 

,00000 

-18.00000 

53 

16 

-7.53220 

-50.39499 

.00000 

-18,00000 

54 

17 

-9.42969 

-50.34838 

,00000 

-18,00000 

55 

la 

- -ia.98323. . 

-50.36485 

,0000.0-. ... 

„ *08. 0.0000 

56 

19 

-12.25515 

-50.36889 

.00000 

-18.00000 

57 

20 

-13.29651 

-50.43786 

.00000 

—18 . 00000 

58 

21 

-14.14911 

-50.49345 

.00000 

-18,00000 

59 

22 

-14*84716 

-50.49759 

*00090 

-18,000.00 

60 

23 

-15.41867 

-50.48889 

.00000 

-18,00000 

61 

24 

-15.88659. 

-50*019  TO 

-ia*nM0jo 

82 

25 

-16.26968 

-50.59025 

,00000 

-18,00000 

63 

26 

r 16. 52334 

-5.31*76836 : 



-18*80000 

64 

i 49 


27 

-16*54014 

-50.76234 

.00000 

-16.00000 

65 

28 

-17.05038 

-50.72418 

.oooco 

-18.00000 

66 

29 

-17.22252 

-50.72116 

.00000 

-18.00000 

67 

30 

-17.36345 

-50.74397 

.00000 

-18.00000 

68 

31 

-17.47884 

-50.73600 

.00000 

-18.00000 

69 

32 

-17.57331 

-50.71086 

.00000 

-18.00000 

70 

33 

-17*65066 

-50.73888 

• Q0Q00 

-18.00000 

71 

34 

-17.71398 

-50.78386 

.00000 

-18.00000 

72 

35 

-17.76583 

-50.82453 

.00000 

-18.00000 

73 

36 

-17.80828 

-50.76659 

..00000 

-18.00000 

74 

37 

17.84303 

-50.82642 

.00000 

-18.00000 

75 

38 

-17.67148 

-50.77776 

.00000 

-16.00000 

76 

39 

-17.89478 

-50.83903 

.00000 

-18.00000 

77 

40 

-17.91385 

-50.75491 

.00000 

-18.00000 

78 

41 

-17.92947 

-50.71252 

.00000 

-18.00000 

79 

42 

-17.94225 

-50.69732 

.00000 

-18.00000 

80 

43 

-17.95272 

-50.75072 

.00000 

-18.00000 

81 

44 

-17.96129 

-50.70645 

.00000 

-18.00000 

82 

45 

-17.96831 

-50.74403 

.00000 

-18.00000 

83 

4b 

-17.97405 

-50.75649 

.00000 

-18.00000 

84 

47 

-17.97876 

-50.76528 

.00000 

-18.00000 

85 

48 

-17.98261 

-50.76587 

.00000 

-18.00000 

86 

49 

-17.98576 

-50.72017 

.00000 

-18.00000 

87 

50 

- *7.98834 

-50.71184 

.00000 

-18.00000 

88 

51 

-17.99045 

-50.66400 

.00000 

-18.00000 

89 

52 

-17.99219 

-50.67188 

.00000 

-18.00000 

90 

53 

-17.99360 

-50.71942 

.00000 

-18.00000 

91 

54 

-17.99476 

-50.76701 

.00000 

-18.00000 

92 

55 

-17.99571 

-50.76317 

.00000 

-18.00000 

93 

5o 

16.98848 

-50.74446 

.01245 

175.00000 

94 

57 

45*63111 

-50.74464 

.07090 

175.00000 

95 

58 

69.08171 

-56.13640 

-10.74591 

175.00000 

96 

59 

88.28144 

-86.31114 

-49.70909 

175.00000 

97 

THIS 

food  is  served 

AT 

TIME  97 

RI  = 

.0000  RO= 

.0000  \IE-  -i 

86.3111 

BAD  SMELL  DETECTED 

IN 

STOR.  2 FOOD 

5 TIME  13 

DURING  SPLIT  NO.  1 

AT 

END  OF  FOOD 

NO.  1 

MAX. 

COUNT  IS  7.4011 

IN  FOOD  NO, 

5 AT  TIME  12 

NSIG=1  NSIGD=1  NTRUzl 


monitor-effectiveness  SUMMARY 
NSz  £ MX-  2 NOC-  El  -Nil-  2 
NSD-  3 NT:  2 ND00=  2 NOllr  2 
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We  make  the  following  comments  about  this  output. 


(i)  The  first  main  section  of  the  output,  ending  with  the  “Intended  Schedule 
of  Splits  and  Blends",  is  fairly  self-explanctory  in  terms  of  the  variable  names  described 
in  Appendix  II.  In  the  “Intended  Schedule  of  Splits  and  Blends"  the  quantifies  captioned 
"Proportion",  are  PR(MS,K),  K = 1,  2,  ...  NJO(MS)  and  MS  is  the  split  number  Also 
the  quantities  captioned  "Handling  Contamination"  and  "Spillage  Weight  Loss"  are  HEVS 
(MS,K)  and  HEWS  (MS,K),  respectively,  for  K = 1,  2,  ...  NJO(MS).  Similarly  for  the 
blends  the  quantities  so-captioned  are  HEVB(MB)  and  HEWB(MB). 

(ii)  In  the  list  of  individual  food  properties  following  the  "Final  Schedule”, 
IE  is  the  time  at  which  the  food  leaves  the  system.  LS  = 0 if  the  food  does  not  end 
in  a split  and,  if  it  does,  LS  is  the  split  number  at  which  the  food  ends.  LB  is  analogously 
defined  for  blends.  The  lines  beginning  "K="  give  the  time,  IACH,  at  which  an  ambient 
temperature  change  occurs,  and  the  new  ambient  temperature,  YAC.  One  such  line  is 
printed  for  each  ambient  temperature  change,  and  no  lines  are  printed  if  the  food  undergoes 
no  change  in  ambient  temperature.  For  example,  food  number  11  suffers  no  ambient 
temperature  change,  while  food  number  13  experiences  two  changes. 

(iii)  In  the  printout  of  the  results  for  the  various  foods,  the  column  headed 
I is  the  number  of  the  time-step,  starting  from  the  time  the  food  entered  the  system, 
while  IAB  gives  the  time  from  the  start  of  the  simulation. 

(iv)  Remember  that  V,  referred  to  as  bacteria  count  in  the  program,  is  really 
th  i natural  logarithm  (i.e.  to  base  e)  of  the  bacterial  density.  Negative  values  of  V merely 
mean  that  the  average  bacterial  density  is  less  than  one  per  gram.  When  V is  large  and 
negative,  it  means  that  the  probability  that  any  bacteria  survives  is  very  low. 

(v)  *«»»mediately  following  the  table  of  temperature  and  bacteria  count  for 
each  food  is  information  abr  it  the  monitoring  systems  and  the  outcomes  of  splits  occurring 
in  the  food.  The  outcome  of  a blend  is  printed  following  the  table  for  the  highest-numbered 
food  entering  the  blend.  If  the  message  "MONREC  SIGNAL"  is  printed,  it  means  that 
the  record-keeping  monitor  system  detected  that  the  food  was  unsafe  (based  on  the  value 
of  TESTI).  If  the  message  "BAD-SMELL  DETECTED"  is  printed,  it  means  that  the 
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domestic-monitor  system  detected  that  the  food  was  unsafe  (based  on  the  values  of  AVI 
and  SD1).  For  example,  in  the  first  trial  of  the  simulation,  a MONREC  SIGNAL  was 
obtained  for  food  number  3 because  RO  = 6 and  TESTI  = 5,  giving  RO  > TESTI.  Also 
in  the  second  trial  the  domestic  monitor  gave  a "BAD-SMELL  DETECTED"  message  at 
blend  number  1,  following  food  number  5. 

(vi)  The  quantities  K and  LIND(K)  give  information  used  in  the  indirectly 
operating  phase  of  the  domestic  monitor  system.  The  quantities  K are  the  numbers  of 
storage  areas  and  LIND(K)  is  an  indicator  that  is  1 if  the  K-th  storage  area  contains  any 
of  the  foods  involved  in  the  split  or  blend  and  0 if  it  does  not.  For  example  in  blend 
number  1,  following  food  number  5,  storage  areas  2,  3 and  4 are  involved  but  storage 
1 is  not.  We  may  verify  this  by  observing  that  foods  2 to  6 are  involved  in  blend  1, 
and  the  Final  Schedule  shows  that  foods  2,  3 and  5 are  in  storage  2,  food  number  4 
is  in  storage  3,  food  number  6 in  storage  4 and  none  of  the  foods  is  in  storage  1. 

(vii)  Following  the  end  of  each  trial  the  maximum  bacteria  count  obtained 
during  that  trial  is  described,  and  the  monitor-system  results  are  summarized.  In  the 
first  trial,  the  record-keeping  monitor  detected  a dangerous  situation  (NSIG=1),  the 
domestic  monitor  system  did  not  (NSIGD=0)  and  the  situation  was  not  really  dangerous 
(NTRU=0).  In  other  words  the  record-keeping  system  gave  a false  alarm. 

(viii)  Also,  the  results  of  the  indirectly-operating  phase  of  the  domestic  monitor 
system  are  described  after  each  trial.  Nothing  is  printed  if  this  monitor  system  detects 
nothing,  but,  if  it  detects  dangerous  situations  (i.e.  a "bad  smell"),  then  messages  are 
printed  out  describing  them.  Such  a message  can  be  observed  at  the  end  of  the  second 
trial  output  above,  where  we  read  that  during  split  1,  at  time  1=13,  a bad  smell  was 
detected  in  food  number  5,  located  in  storage  2.  Referring  back  to  food  number  5, 
we  see  that  V = 7.34  at  I = 13,  and  since  AV2  = 7.20,  SD2  = 0,  the  count  was  high 
enough  to  exceed  the  threshold  level,  given  in  this  case  by  AV2.  In  the  first  trial,  the 
corresponding  value  of  V is  7.075, which  was  less  thanAV2,  and  therefore  no  bad  smell 
was  detected. 

(ix)  The  monitor-effectiveness  summary,  given  at  the  very  end  of  the  print-out, 
shows  that  the  record-keeping  monitor  gave  detections  in  all  of  the  5 trials  (NS=5)  that 
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were  run,  but  only  in  2 cases  was  there  truly  a dangerous  situation  (NT=2),  hence  that 
monitor  system  gave  false-alarms  in  3 of  the  5 cases  tried.  The  domestic  monitor  detected 
danger  in  3 out  of  5 cases  (NSD=3),  correctly  classifying  2 cases  where  there  was  no 
danger  (NDOO=2)  and  two  cases  where  there  was  danger  <N D 1 1=2) , but  giving  one  false 
alarm. 
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Appendix  IV:  Program  Listing 


The  main  part  of  this  Appendix  is  a listing  of  the  Fortran  V source  program,  called 
BACTRAC,  which  performs  the  simulation.  There  is  quite  a lot  of  explanation  of  the 
program  in  this  listing  by  means  of  comment  lines,  i.e.  lines  of  text  beginning  with  a 
C in  column  1.  Since  it  is  moderately  self-explanatory,  we  give  only  the  following  brief 
comments  as  guidance: 

(i)  The  heart  of  the  program  is  the  big  loop  that  begins  at  line  206  and  ends 
at  line  464.  Each  passage  through  this  loop  does  one  trial  of  the  entire  simulation. 
Everything  that  precedes  this  big  loop  is  reading  in  of  the  data  and  initialization  of  some 
parameters.  Within  this  loop  is  a smaller  (but  still  large)  loop,  extending  from  line  291 
to  line  420.  Each  passage  through  this  loop  does  the  calculations  for  one  food,  primarily 
by  calling  the  subroutine  BACGRO  at  line  302. 

(ii)  The  subroutine  SPLEND  does  the  calculations  involved  in  the  split  or  blend 
that  occurs  after  foods  have  been  in  storage.  It  is  called  twice  in  the  program,  at  lines 
351  and  390.  The  number  IOP  determines  whether  SPLEND  does  a split  (IOP=1)  or 
a blend  (IOP=2)  calculation. 

(iii)  The  random  numbers  used  in  the  program  are  generated  by  the  system 
subroutine  RANDN,  which  is  called  at  three  different  places  in  the  program,  lines  210, 
485  and  543.  RANDN(X,  N,  A,  S)  generates  N random  numbers,  normally  distributed 
with  mean  A and  variance  S2,  and  stores  them  in  the  array  named  X.  If  the  program 
is  run  on  a system  that  does  not  have  a random  number  generator  of  this  type  in  the 
system  subroutines,  the  user  must  write  one  and  add  it  to  the  subroutines  of  the  program. 
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1* 

2* 

3* 

4 ♦ 

b* 

b* 

7* 

b* 

9* 

1U* 

11* 

lc* 

13* 

14* 

lb* 

lb* 

17* 

la* 

19* 

2u* 

21* 

2<r* 

23* 

24* 

2b* 


C BACTRAC  PROGRAM 


THE  INSTRUCTIONS  OP  THE  RECORD-KEEPING  MONITOR  SYSTEM  ARE  MARKED 
BY  MONREC  IN  COLUMNS  73-78.  THE  INSTRUCTIONS  OF  THE  DOMESTIC 
MONITOR  SYSTEM  ARE  MARKED  BY  MONDOM  IN  COLUMNS  73-78. 

THE  INDEX  J ALWAYS  REFERS  TO  THE  VARIOUS  FOODS. 


DIMENSION  X ( 140 ) *IC(20) > I E ( 20 ) > W ( 20 ) » Y (20.250  > * V (20.250 ) » I CD (20) . 

1 WD(20) .YD (20) *VD(20)  . ISI20)*  10(20).  JSI (20) . JBO(20) . ISD 

2(20) . IBD (20 ) > JSO (20 . 20 ) . JBI ( 20 . 20 ) . TITLE ( 12 ) .LEGEND ( 12.20 ) 

DIMENSION  NJO(20) *NJI (20) »LS(20) .LB (20) , Y A ( 20 ) « P ( 20 . 6 ) »Q(6) 
DIMENSION  I AC (20.5) .IACD120.5) »YAC(20.5) .NATC (20) .LIND (20 ) 
DIMENSION  IO(20)»YO(20).YE(20).VE(20)» ISER(20) .MST (20) . YST (20) 
DIMENSION  JJ(5) »YY (5) »VV(5) • MW ( 5 ) »PR(5»5) »PRP(5) »CONT(20) »YSTD(20) 
DIMENSION  HE VS (2 0.5) »HEWS(20«5> . HEVB ( 20 ) .HEWB (20 ) »HVS(5) »HWS(5) 
DIMENSION  IO(20).Yl(20)»VI(20)»Wl(20)»lND(20).lACH(20»5) 

DIMENSION  IREF(20) .VTRAP<20) .VREC (20.250) .JREC (20) .IV(20) 

DIMENSION  Hi (20) .R0(20) ,LRS(20) .LRB120) ,JQ(20) . I VM ( 20 ) . VMAX ( 20 ) 

200  FORMAT  (3H0J=»I2»2X.12A6) 

201  FORMAT  ( 5X , 4 ( 12 . 2X ) . 2X . 3 ( E 1 2 . 5 . IX ) , 2X . El2 . 5 ) 

202  FORMAT (1H0.5X.11HFOOD  NUMBER . 1 3 , 1 X 12 A6 ) 

203  FORMAT  (1H0»5X. ’FINAL  SCHEDULE*) 
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<26* 
27* 
28* 
29* 
30* 
31* 
32* 
33* 
34* 
3b* 
36* 
37* 
38* 
39* 
40* 
*1* 
42* 
43* 
44  * 
4b* 
4o* 
47* 
48* 
49* 
bO* 
bl* 
b2* 
b3* 

b4* 

bb* 

b6* 

b7* 

b6* 

b9* 

bQ* 

bl* 

o2* 

63* 

64* 

bb* 

66* 

b7* 

b8* 

b9* 

70* 

71* 

72* 

73* 

74* 

7b* 

76* 

77* 

76* 

79* 

60* 

dl* 

82* 


204  FORMAT ( • FOOD' , 13, ' STARTS  AJ  ls',I3»»  ANO  ENDS  AT  SPLIT'. I3»'  STo 
1R.  N0.=«#I2> 

205  FORMAT ( • FOOD', 13,'  STARTS  AT  Is', 13,'  AND  ENDS  AT  BLEND' 1 13. ' STO 
1R.  N0.=**I2) 

206  FORMAT  (•  FOOD', 13.'  STARTS  AT  1 = * • 13* ' AND  CONTINUES  TILL  THE  END 
l»Is',l3»'  STOR.  N0.=*#I2) 


207  FORMAT  (lHQ.SX, 'INTENDED  SCHEDULE  OF  SPLITS  AND  BLENDS') 

208  FORMAT  POSPLIT  NUMBER', 12, » IN  FOOD  NUMBER'»I3»'  OCCURS  AT  TIME'. 
1 14  » ' WITH  STD  DEV  ' .13) 

209  FORMAT  (•  THE  NUMBERS  AND  PROPORTIONS  OF  THE  RESULTING  FOODS  ARE') 

210  FORMAT  ( 3X . 'FOOD  NUMBER'. I3»*  PROPORTION' ,F8. 5) 

211  FORMATdHO, 'BLEND* . 12. • PRODUCES  FOOD'. 13.'  AT  TIME'»I4»'  WITH  S.D 
1.  ' » 13. • FROM  FOODS' ) 

212  FORMAT  (IX, 17(13. IX) ) 

213  FORMAT  (4X,'I'.6X»'Y=TEMP',5X.'V=BAC.  COUNT' ,7X. ' A* ,6X, ' YAsAMBIENT 
1 TEMP' ,)X, 'IAB« ) 

214  FORMAT  CoTHlS  FOOD  ENDS  IN  SPLIT*. 13.'  AT  TIME.  I='.l4.».  INPUT 
1 TO  SPLIT  HAS') 

215  FORMAT  (3X,*TEMP=Y=',F8.3.3X'BAC.  COUNTsVs* ,F8.3. ' WEIGHTsWs* ,F8.3 
1> 

216  FORMAT  (•  OUTPUT  FROM  SPLIT  IS  AS  FOLLOWS') 

217  FORMAT  (3X*FOOO  NO.  '.12.'.  TIMES', I3»'»  Y=',F8.3.'»  Vs,,F8.3,'»  W 


1=' .F8.3) 

216  FORMAT  COTHIS  FOOD  ENDS  IN  BLEND  NO.*.I3.»,  INPUT  TO  BLEND  HAS') 
219  FORMAT  (3X, 'FOOD  No.  '.12,',  TIME=».I3.».  r='.F8.3.'.  Vs',F8.3,'. 
1WT.='.F8.3) 

OUTPUT  FROM  BLEND  HAS') 

(•FOOD  NO. ='  I2»'.  TIMES'. 13. '.  Ys',F8.3»'»  Vs',F8.3.*,  W 


THIS  FOOD  IS  SERVED  AT  TIME'. 13) 

This  FOOD  IS  STORED  UNTIL  THE  FINAL  TIME.'. 14) 

FOOD', 13.'  STARTS  AT  Is', 13.'  AND  IS  SERVED  AT  Is', 13.' 
s*. 12) 

X,10(F10.5. IX)  ) 

NJs',12.'  NMSs',12,' 

9.6,'  NREPs ' » 14  ) 
lCs',13.'  Ws* ,F9.4» ' 

.s'  .12) 

ICD='.I2,'  WDs» ,F8.4, • 


220 

FORMAT 

( • 

22 1 

FORMAT 

■3 

lT.s* ,F0.3 

222 

format 

(1 

223 

format 

( • 

224 

format 

. 

225 

FORMAT 

( • 

1STOR.  NO. 

22b 

format 

(1 

227 

format 

( • 

1'  SDVs • 

, F 

226 

format 

( • 

1 STOR. 

NO 

229 

FORMAT 

( ' 

1.'  CONT 

- f 

230 

FORMAT 

( • 

231 

format 

( • 

232 

format 

( ’ 

233 

format 

( • 

234 

format 

( ' 

235 

FORMAT 

1) 

( • 

236 

format 

( ' 

237 

format 

( • 

238 

format 

( • 

239 

format 

( ' 

240 

format 

( • 

241 

format 

( • 

242 

format 

(• 

243 

format 

( • 

244 

format 

( • 

NMBs',12,'  X ( 1 ) s ' , F7 . 0 » * lMAXs',14, 
Yls* ,F8.4» * ;is',F7.4,»  ISERs • , 1 3 » • 
YDs ' , F8 • 4 , • VDs» ,F8.4, ' NATCs',l2 


P(J»K) » ,b(lX,F9.5)  ) 

K=' »I2. • IACs.,13,.  IACD: 

KS • » I2» • IACHs',13,’  YAC: 
lCs',13.'  IE=' » 13. * LSs' » I2» ' 

*s',F9.4.’  Yls' ,F9.4. ' VIs'.F9.4,'  YAs',F9.4) 
ISTORAGE  UNIT', 13.'  IS  AT  TEMP', Fa. 3.'  WITH  STD 


',12,'  YAC=' »F9.4) 
• » F9 . 5 ) 


Lbs » , 12 . ' NATOs* ,12) 

DEV  ,F8.4 


bad  SMELL  DETECTED') 

BAD  SMELL  DETECTED  JUST  BEFORE  SERVING') 

AVIs', F8. 3.'  SDls»,F8.3»»  AV2s',f8.3»*  SD2s',F8.3) 

BAD  SMELL  DETECTED  IN  STOR. ',12.'  FOOD*. 12.'  TIME*. 13) 
NSIG=  * , 1 1 , ' NSIGDS'.II.’  NTRUs',H) 

NSs ' , 13. ' NTs', 13, » NOOs',13,'  Nils', 13) 

DURING  SPLIT  NO. »,I3»'  AT  END  OF  FOOD  NO.', 13) 

DURING  BLEND  NO.»,I3»*  AT  END  OF  FOOD  NO. '.13) 

MAX.  COUNT  IS' »F9.4» ' IN  FOOD  NO. ',13,'  AT  TIME', 13) 
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* 


83* 

245 

FORMAT 

<• 

NSD='.I3.'  NTs*.  13,*  ND00=*»l3.*  NDU='.I3) 

84* 

246 

FORMAT 

( * 

HANDLING  CONTAMINATION* .F10. 5. * . SPILLAGE  WEIGHT 

LOSS 

1 

85* 

1.F10.5) 

86* 

247 

format 

<* 

MONREC  SIGNAL*) 

87* 

248 

FORMAT 

(’ 

MRN0=*.I3.'  IPTRD=*»U.*  NST0R=»,I2) 

88* 

249 

FORMAT 

( * 

MONITOR-EFFECTIVENESS  SUMMARY*) 

89* 

300 

FORMAT 

< * 

Ql=' *F8.4»  * Q2=»*Fa.‘i.  • Q3=*  *F8,4.  * TEST1=*  .F8.4.  * 

TRT= 

f 

90* 

1.F8.4) 

91* 

301 

FORMAT 

(• 

RI='.F10.4) 

MONREC 

92* 

302 

FORMAT 

<* 

RI=*.F10.4.*  R0=».F10,4.*  VE=*.Fl0.4  ) 

MONREC 

93* 

100 

FORMAT 

(12A6) 

94* 

101 

FORMAT 

() 

95*  C *******************************  *****#********+**+***************************** 

96*  C 

97*  C WE  READ  IN  PROPERTIES  OF  THE  FOOD  PACKETS 

98*  C 


99* 

100* 

101* 

102* 

103* 

104* 

10b* 

106* 

107* 

108* 

109* 

110* 

111* 

112* 

1 113* 
114* 
115* 
116* 
117* 
118* 
119* 
120* 
121* 
122* 
123* 
124* 
125* 
126* 
127* 
128* 
129* 
130* 
131* 
132* 
133* 
134* 
13b* 
136* 
157* 
138* 
139* 


C*** *********************************************************** **************** 
WRITE  (6*222) 

READ  (5.100)  (TITLE  (K)rK=lil2) 

WRITE(6.100)  (TITLE  (K).Ksl,12) 

READ  (5.101)  NJ.NMS.NMB.XU) , IMAX .SDV .NREP ,NRNO , IPTRD.NSTOR 

WRITE  (6.227)  NJ.NMS.NMB.X(l) . IMAX.SDV.NREP 

WRITE  (6.248)  NRNO. IPTRD.NSTOR 

READ  (5.101) AVI. SD1.AV2.SD2 

WRITE (6. 238) AV1.SD1.AV2.S02 

READ  (5.101)  Q1 .02.03  .TEST1.TRT  MONREC 

WRITE(6.30Q>  Q1.Q2.Q3  .TEST1.TRT  MONREC 

DO  10  Jrl.NJ 

READ  (5.100)  (LEGEND  (K. J) ,K=1. 12) 

WRITE  (6.200)  J. (lEGENO(K, J) .K=1.12) 

READ  (5.101)  IC(J),W(J),Y(J.1),V(J.1).ISER(J).MST(J) 

WRITE  (6.228)  IC ( J ) . W ( J) . Y ( J , 1 ) , V ( J. 1 ) . ISER ( J ) .MST (0 ) 

READ  (5.101)  ICD(J).  WD ( J) . YD ( J ) . VD ( J ) , NATC ( J) ,CONT ( J) 

WRITE (6 . 229 ) ICD ( J) » WD ( J) , YD ( J ) . VD ( J) .NATC (0) . CONT ( J i 

READ  (5.101)  <P(J.K).K=1.6> 

WRITE (6. 230) (P( J.K) .K=l ,6) 

READ  (5.101)  RI(J) 

WI (J)=W(J) 

YI(J)=Y(J.l) 

VI(J)=V(J. 1) 

IQ(J)=IC(J) 

N=NATC(J) 

IF  (N.LT.l)  N=1 

READ  (5.10I.ERR=10)  ( I AC ( J . K ) . I ACD ( J. K ) , YAC ( J *K ) ,K=1 .N ) 

DO  46  K=1.N 

46  WRITE  (6.231)  K, I AC ( J. K) . I ACD ( J.K ) . YAC ( J .K ) 

IF  (IC(J).EQ.l)  WRITE  (6.301)  RI(J)  MONREC 

10  CONTINUE 

READ  (5.101)  (YST(KS) .YSTD(KS) .KS=1.NST0R) 

DO  47  KS=1,NST0R 

47  WRITE  (6.235)  KS.YST(KS) .YSTD(KS) 

WRITE  (6.207) 

C******** ********************************************************************** 

c 

C WE  READ  IN  THE  PROPERTIES  OF  THE  SPLITS 

C 

C **************************************************************** ************** 
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140* 

l4l* 

142* 

143* 

144* 

145* 

146* 

147* 

148* 

149* 

150* 

151* 

152* 

153* 

164* 

155* 

15b* 

157* 

156* 

159* 

160* 

lbl* 

162* 

163* 

164* 

165* 

166* 

167* 

166* 

169* 

170* 

171* 

172* 

173* 

174* 

175* 

176* 

177* 

178* 

179* 

180* 

181* 

162* 

183* 

184* 

185* 

186* 

187* 

188* 

189* 

190* 

191* 

192* 

193* 

194* 

195* 

196* 


C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

c 


DO  38  J-liNJ 


38  IND(J)=0 

00  11  MS-1 *NMS 
READ  (5.101)  NJO(MS) 

N=NJO(MS) 

READ  (5.101)  IS(MS)  » JSI  (MS)  . ISO(MS) . (vJSOIMS 
REAO  (5.101)  (HEVS(MS.K) .HEWS(MS.K)  .Krl.N) 
DO  28  K=1.N 
JOY-JSO(MS.K) 

28  IN0(J0Y)=1 

WRITE  (6.208)  MS. JSI (MS) « IS(MS) . ISD(MS) 
WRITE  (6.209) 

WRITE  (6.210)  (JSO(MS.K) .PR(MS.K) .K=1.N) 
WRITE (6. 246)  (MEVS(MS.K) .HEWS(MS.K) .K=1,N) 
11  CONTINUE 


.K) .PR(MS. K) ,K=1,N) 


we  read  in  the  properties  of  the  blends 


DO  12  MB=1»NMB 
READ  (5. lui < NJI(MB) 
N=NJI (MB) 


READ  (5.101)  IB(MB).JBO(MB).IBD(MB)»(JBI(MB.K).K= 
READ  (5.101)  HEVB(MB) .HEWB(MB) 

JOY=JBO(MB) 

IND ( JOY ) - 1 


l.N 


) 


WRITE  (6.211)  MB. JBO(MB) . IB(MB) . IBD(MB) 
WRITE  (6.212)  (JBI(MB.K),K=1,N) 

WR I TE ( 6 » 246 ) HEVB(MB) .HEWB(MB) 

12  CONTINUE 


we  determine  for  each  food  the  quantities  lsijj  an,?  lbjj),  which 
specify  whether  the  j-th  food  ends  in  a split  or  a blend. 


00  13  Jrl.NJ 

LS(J)=0 

LB(J)=0 

DO  14  MS=1,NMS 

NKO-N.'O  (MS) 

DO  15  Krl.NKO 

15  IF  (JSI (MS)  .EQ.J)  LS(J)=MS 
14  CONTINUE 

DO  16  MB=1.NMB 
NKI=NJI (MB) 

DO  17  Krl.NKI 

17  IF  ( JBI ( MB. K) .EQ.J)  LB(J)=MB 

16  CONTINUE 
13  CONTINUE 


C********  **♦*******♦♦***♦*♦************************, **#**********<^^,m^^^ 

C HERE  WE  INTRODUCE  RANDOMNESS  IN  INITIAL  VALUES  AND  IN  SCHEDULING. 

C ALSO  THE  BIG  LOOP  THAT  ENDS  AT  INSTRUCTION  NO.  50  DOES  THE  REPETI- 

C TIONS  OF  THE  ENTIRE  SIMULATION. 

C 
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197* 

C**«* 

lsa* 

N0Q  = 0 

199* 

N11  = 0 

200* 

NS=0 

2ul* 

NSD=0 

202* 

ND00  = 0 

203* 

NDllsO 

2 04  * 

NT  = 0 

2ub* 

NREF=NMS*NM6 

2ub* 

DO  50  M=1.nREP 

20  7* 

NSIGD=0 

2uo* 

NSIG=0 

209* 

NTRUrO 

clG* 

CALL  RANDN  ( X . NRNO » 0 . r 1 . ) 

^11* 

DO  29  J=l.NO 

i.'12* 

IF  (NATC(J) .LT.11  GO  TO  43 

2l  3* 

N=NATC ( J) 

2W* 

NNR  = NMS*-NMb  + 4*NJ*l 

210* 

DO  44  K=1*N 

2 lo* 

IACH( J,K)=IAC(J.K)+IACDC J*K)*X(NNR) 

217* 

NNR=NNR*1 

216* 

44 

CONTINUE 

219* 

43 

IF  ( IC(J)  .(jT.l)  GO  TO  42 

220* 

GO  TO  39 

2£1* 

42 

IF  ( IND(J) .EO.l)  GO  TO  29 

222* 

ICl J)=IQ( J)+ICD(J)*X(J*3*NJ) 

223* 

39 

IA=IC(J) 

224* 

W(u)=WI(J)+WD(J)*X(J) 

22b* 

Yl Jr IA)=YI (J)+YD(J)*X( J*NJ) 

22  b* 

V<  Jr  IA)=VI(J)*VO(J)*X(J42*NwI) 

227* 

29 

CONTINUE 

226* 

DO  32  MS=1,NMS 

229* 

N=NJO(MS) 

2 JO* 

DO  33  K=lrN 

2 J 1 * 

J=JSO(MSrK) 

2 Jc« 

IC( J)=IQ( J)*ISD(MS)*X(MS*4*NJ) 

2 Jj* 

33 

CONTINUE 

2 J4  * 

IS(MS)=IC ( j) 

2 Jb* 

IREF(MS)~IS<MS) 

2 Jb* 

32 

CONTINUE 

2 J7* 

DO  37  MB=I,NME 

2 JU* 

J=JB0(MB> 

-39* 

IU J)=IQ(J)+I3D(MB)*X<NMS+4*NJ*MB) 

240* 

IB(MB)=IC(j) 

241* 

IREF (NMS*Mb)=IB(MB) 

37 

CONTINUE 

24  J* 

WRITE  (6*203) 

244* 

DO  34  J—  1 1 

240* 

IF (LS< J) .EQ.O)  GO  TO  35 

LHtl* 

WRITE  <b#204;  JrIC(J) rL5(0)rM5T(J) 

2 4 7 ♦ 

GO  TO  34 

246* 

3b 

IF  (LB ( J) .EQ.OJ  GO  TO  36 

C49* 

WRITE  (6*205)  JrIC(J) rLB<J) rMST(J) 

2b  0* 

GO  TO  34 

2hl* 

36 

IF  (ISER(J) .GE.IMAX)  WRITE  <6*2061  Jr IC  ( J) , IMAX rMST ( J ) 

202* 

IF  ( ISERt J) .LT.IMAX)  WRITE  (6r225>  Jr IC ( J) , ISER ( J) rMST < J > 

2b  J* 

34 

CONTINUE 

MONREC 

MONREC 

MONREC 

MONDOM 

MONDOM 

MONDOM 

MONDJM 

MONDOM 

MONREC 


MONDOm 


MONDOM 
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254* 

C* 

••*1 

255* 

c 

256* 

c 

HERE  WE  CALCULATE  THE  TIME-STEP  AT  WHICH  EACH  FOOD  TERMINATES 

257* 

c 

256* 

c* 

• ** 

259* 

00  IB  Jsl.NJ 

260* 

IF  (LS(J).LT.l)  00  TO  19 

2ul* 

MS=LS(J) 

262* 

!E(J)=IS<MS> 

2b3* 

GO  TO  21 

269* 

19 

IF  (LB(J).lT.I)  60  TO  20 

2o5* 

M6=LB( J) 

266* 

IE(J)=IB(MB> 

2b7* 

GO  TO  21 

2b8* 

20 

IE(0)=IMAX 

269* 

IF  (ISER(J) .LT.IMAX)  IE(J)=ISER(J) 

270* 

21 

WRITE  (6*200)  J*  (LEGEND (K ( j ) , K = 1 * 12) 

271* 

IA=IC( J) 

272* 

MPsMST ( J) 

27  3* 

YA(J)=YST(mP)+X(NNR)*YSTD(MP) 

274* 

NNR=NNR*1 

275* 

WRITE  (6*233)  IC  ( J) * IE ( J) *LS ( J) *LB ( J) .NATC l J) 

276* 

WRITE  (6.234)  W(J) *Y(J.IA)  ,V<J,IA) ,YA(J) 

277* 

N=NATC( J) 

278* 

IF  (N.LT.l)  GO  TO  18 

279* 

DO  45  K=1.N 

280* 

45 

WRITE  (6*232)K.  IACH(J.K) *YAC(J*K) 

281* 

18 

CONTINUE 

282* 

c* 

**♦ 

******* 

2o.J* 

c 

284* 

c 

HERE  BEGINS  THE  MAIN  LOOP. WHICH  CALCULATES  TEMPERATURE  AND  BACTER- 

265* 

c 

IAL  GROWTH  IN  EACH  FOOO  AND  THEN  MAKES  THE  CALCULATIONS  ARISING 

28b* 

c 

OUT  OF  THE  SPLIT  OR  BLEND  AT  WHICH  THAT  FOOD  PACKET  TERMINATES. 

287* 

c 

266* 

c* 

*** 

289* 

NNR=NMS*NMB*4*NJ*2 

290* 

L = 0 

MONDOM 

291* 

DO  52  J=1*nU 

292* 

II=1C(J) 

295* 

DO  23  K = 1 ,6 

294* 

23 

Q(K)=P(J,K) 

295* 

N=NATC(J) 

296* 

DO  24  K=1.N 

297* 

YO(K)=YAC(j.K) 

298* 

24 

IO(K)=IACH(J.K) 

299* 

WRITE  (6*202)  J. (LEGEND  (K.J) *K  = 1 * 12) 

3UU* 

IF  (IPTRD.NE.l)  WRITE  (6*213) 

301* 

YAA=YA( J) 

302* 

41 

CALL  BACGRO  (Q.II.IE(J) .lO.YA(J) *Y(J.II ) . YO.YE(J) * V ( J* 1 1 ) .SDV.NAT 

303* 

1C(U) ,X(NNR) *VE(J) *IREF.NREF*VTRAP*VMAX(J) , IVM(J) *CONT(J) ) 

304* 

YA(U)=YAA 

305* 

NNR=NNR*1 

30b* 

00  54  Krl.NREF 

MONDOM 

307* 

IZ=IREF(K) 

MONDOM 

308* 

VREC ( J* IZ)=VTRAP(K) 

MONDOM 

309* 

54 

CONTINUE 

MONDOM 

310* 

IF  (ISER(J) .GE.IMAX)  GO  TO  48 

60 


Ml* 

312* 

313* 

314* 

31b* 

31b* 

317* 

316* 

319* 

320* 

321* 

322* 

323* 

324* 

325* 

32b* 

3*7* 

320* 

329* 

330* 

331* 

332* 

333* 

334* 

335* 

33b» 

337* 

330* 

339* 

340* 

341* 

342* 

343* 

344* 

34b* 

34b* 

347* 

34b* 

349* 

3b0« 

351* 

352* 

353* 

3 J4* 

355* 
356* 
357* 
350* 
3d9« 
3o0  * 
3d  1 * 
3q2  * 
3o3* 
304* 
3o5* 
3do* 
3b  7* 


C 

C 

c 

c 


c 

c 

c 

c 


WRITE  (6*223)  IE  < J ) 

CALL  SNIFF  ( VE  ( J )**  (J )* AVI , SOI  * X ( NNR )* NOSE ) MONDOm 

IF  (NOSE.EQ.l)  WRITE  (6*237)  MONDOm 

nsigd=nsigd*nose 

NNR=NNR+1 
GO  TO  49 

46  IF l (LS( J)+lB< J) ) .EQ.O)  WRITE  (6.224)  IMAX 
49  IT=IE(J) 

Y(J. IT)rYE(J) 

V( J*IT)=VE(J) 

YCOMP=YA( J) 

IF  (N.EQ.O)  GO  TO  57 
DO  5b  K = 1*N 

56  YCOMPzAMAXl (YCOMP. YO(K) ) 

57  RO(J)=0. 

IF  ( YCOMP. lT.QI)  RO(J)zRI(j) 

IF  ( ( YCOMP. GE.Q1 ) .AND. (YCOMP. LE.Q2) ) R0< J)=RI ( J ) +Q3* ( YC0MP-Q1 ) * ( I T 
1- I I ) MONREC 

WRITE  (6*302)  RI  ( J) » R 0 ( J ) *VE( J)  MONREC 

IF  ( RO ( J ) , GE . TEST  1 ) NSIG=NSIG+1  MONREC 

IF  ( RO ( J) .bE.TESTl)  WRITE  (6*247) 

From  here  to  25  DESCRIBES  WHAT  HAPPENS  TO  A FOOD  THAT  ENDS  IN  A 
SPLIT. 

IF  l LS ( J) .tO.O)  GO  TO  25 
Mb=LS(J) 

NzNJO(MS) 

CALL  SNIFF  ( VE ( J ).  w ( J >. AV 1 , SD1 . X < NNR )* NOSE ) MONDOm 

NSIGD=NSIGD+N0SE 

NNR=NNR*1 

I0P=1 

WRITE  (6*214)  MS, IT 

WRITE  <6.2l5.ERR=40>  YE ( J ) , VE ( J ) * W ( J ) 

IF  ( NOSE • EG • 1 ) WRITE  (6.236)  MONDOm 

40  DO  26  K=1,N 

PRP(K)=PR(mS*K) 

HVS(K)=HEVS(MS.K) 

HwS(K)=HEWS(MS*K) 

2b  JJIK)ZJSO(MS.K) 

call  SPLEND  ( I OP* JSI (MS) *N,YE< J) *VE(J>  *W(J) . JJ. YY * VV.WW.PRP.X (NNR ) 

1 .HVS.HwS.NSTOR.MST.LIND) 

NNR=NNR+1 
WRITE  (6.216) 

DO  27  K=1,N 
JO=JJ(K ) 

RI(JO)=RO(J)  MONREC 

1 1 = I C ( JO) 

Y(J0.II)=YY(K) 

V ( JO  , 1 1 ) =w  (K  > 

W ( JO ) =WW ( K ) 

27  WRITE  (6.217)  JO  * 1 1 * Y ( JO  * 1 1 ) . V ( JO . 1 1 ) * W ( JO ) 

GO  TO  22 

FROM  HERE  TO  22  DESCRIBES  WHAT  HAPPENS  WHEN  SEVERAL  FOODS  END  IN 
A BLEND  SIMULTANEOUSLY. 
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3b8* 

2b 

IF  (LB(J).EO.O)  GO  TO  52 

3b9* 

MB=LB( J) 

370* 

N=NJI (MB) 

371* 

JMAX=0 

372* 

00  30  K=1*N 

373* 

30 

JMAXrMAXO  ( JMAX, JBI (MB.K) ) 

37*+* 

IF  (J.NE.JMAX)  GO  TO  52 

37b* 

I0P=2 

37b* 

WHITE  (6*218)  MB 

377* 

JU=JBO(MB) 

378* 

B=0. 

MONREC 

379* 

DO  31  K=1  *N 

300* 

JJ(K)=JBKmB»K) 

3ol* 

JR=JJ(K) 

3b2* 

B=b+RO( JR) 

MONREC 

3o3* 

IT=IE( JR) 

3o*+* 

YY(K)=Y(JR,IT) 

385* 

VV(K)=V( JR, IT) 

30o* 

mm (K ) =M ( JR ) 

3o7* 

WRITE  (6*219)  JR.IT.YY(K) *VV(K) , WW ( K ) 

38b* 

31 

CONTINUE 

389* 

RI ( JU)=B/N 

MONREC 

390* 

CALL  SPLENO  ( I OP . JU » N . YU  * VU . WU . JJ * YY * VV , WW , PRP * X (NNR ) *HEVB (MB ) .HEW 

391* 

IB (MB ) * NSTOR  * MST  * LIND ) 

392* 

WHITE  (6*220) 

393* 

NNR=NNR*1 

39*+* 

I 1 = I C ( JU ) 

395* 

Y(JU.II)=YU 

396* 

V(JU,II)=Vo 

397* 

wuu;;wu 

398* 

WRITE (6 *221)  JU*I1 *Y(JU,II) ,V(JU*II) *W(JU) 

399* 

CALL  SNIFF  (VU*WU*AV1»S01*X(NNR) .NOSE) 

MONDOM 

*400* 

NNR=NNR*1 

*401* 

NSlGO=NSIGb*NOSE 

4u2* 

IF  (NOSE.EQ.l)  WRITE  (6.236) 

MONOOM 

403* 

22 

CONTINUE 

404* 

C 

405* 

C 

the  loop  ending  at  si  determines  the  foods  checked 

BY  THE  DOMESTIC 

40b* 

C 

MONITOR  SYSTEM  DURING  THE  SPLIT  OR  BLEND  BEING  DONE 

. • 

4U7* 

c 

408* 

DO  51  JIMr i , NJ 

40  )* 

K=MST( JIM) 

41  J* 

IF  (LIND (K ) • EG , 0 ) GO  TO  51 

411* 

IF  ((II.LE.IC(JIM)).OR.(II.GE.IE(vJIM))>  GO  TO  51 

412* 

L=L  + 1 

413* 

I V (L ) = 1 1 

414* 

JREC(L)=JIm 

415* 

LHS(L)=LS( j) 

4 lo* 

LRB(L)=LB(J) 

417* 

JO(L)=J 

418* 

51 

CONTINUE 

419* 

NL=L 

MONDOM 

420* 

52 

CONTINUE 

42l* 

c 

422* 

c 

the  loop  ending  at  53  is  also  part  of  the  domestic 

MONITOR  SYSTEM. 

423* 

c 

IT  CHECKS  ALL  FOODS  INVOLVED  IN  THE  SPLIT  OR  BLEND 

BEING  DONE*  TO 

424* 

c 

SEE  IF  THE  COOKS  SMELL  ANYTHING  BAD. 

62 


425* 

<♦26* 

427* 

428* 

429* 

*430* 

431* 

*432* 

433* 

434* 

435* 

436* 

437* 

438* 

439* 

440* 

441* 

4 4 2 * 

443* 

444* 

445* 

440* 

447* 

44  6* 

449* 

430* 

431* 

452* 

453* 

454* 
435* 
450* 
437* 
438* 
459* 
460* 
401* 
462* 
463* 
4b4* 
4q5* 
4oO  * 
4q7* 
466* 
4b9* 
470* 
471* 
47^* 
473* 

4 74* 
473* 

4 7 o * 

4 7 7* 
476* 
479* 
460* 
461* 


C 

00  53  L=1*NL 
I=IV(L) 

J=JR£C(L) 

VSM=VREC(J,I) 

CALL  SNIFF  < VSM, w ( J ) * AV2 > SD2 * X (NNR ) , NOSE ) 

NNR=NNR+1 

IF  (NOSE.NE.l)  GO  TO  53 
WRITE  (6*239)  MST ( J ) * J « I 

IF  (LRS(L).GT.O)  WRITE  (6*242)  LRS(L)*JQ(L) 

IF  (LRB(L) .GT.O)  WRITE  (6*243)  LRB (L ) * JQ (L ) 

NS1G0=NSIG0*1 
53  CONTINUE 
C 

C FROM  HERE  TO  50  WE  CALCULATE  THE  MAX  BACTERIA  COUNT  AND 
C INFORMATION  FOR  THE  MONITOR-EFFECTIVENESS  SUMMARY. 

C 

VMD  = 0 . 

DO  55  J=1.NJ 

IF  ( VMC.GE.VMAX( J) ) GO  TO  55 
VMD=VMAX ( J) 

JM  = J 

IMM=IVM< J) 

55  CONTINUE 

IF  (VMD.GE.TRT)  NTRU=NTRU+1 
WRITE  (6*244)  VMD.JM.IMM 
X ( 1 ) =X (NNR ) *10000000 
IF  (NSIG.GE.l)  NSIG=1 
IF  (NTRU.GE.l)  NTRU=1 
IF  (NSIGD.GE.I)  NSIGD=1 
WRITE  (6*240)  NS IG . NSIGD  * NTRu 
WRITE  (6*222) 

NS-NS*NSIG 
NSD-NSO+NSIGD 
N r zNT  + NTRU 

IF  ( (NSIG.EQ.O) .AND. (NTRU. EQ.O)  ) NOOzNOO+l 
IF  l (NSIG.EQ.l)  .AND.  (NTRU.EQ.l)  ) N11=N11+1 
IF  ( (NSIGD. EQ.O)  .AND.  (NTRU. EO.O) ) NDOO=NDOo*l 
IF  ( (NSIGD. EG. 1) .AND. (NTRU.EQ.l)  ) ND11=ND11  + 1 
50  CONTINUE 

WRITE  (6*249) 

WRITE  (6*241)  NS  *NT  * NOO  *NJ  1 
WRITE  (6*245)  NSD*NT.ND00,ND11 
STOP 

C*********************************************************#**f*#****$*$*$rt*$$^$ 

C THIS  IS  THE  END  OF  THE  MAIN  PROGRAM*  WE  NOw  START  THE  LISTING  OF 

C THE  SUBROUTINES. 

c the  subroutine  bacgro  does  the  calculation  of  temperature  change 

C AND  BACTERIAL  GROWTH  FOR  EACH  FOOD. 

c 

C******«****t*»««**«M******«****M******»****«#*******,*»t***M*,M4,,t,,,M(t 

SUBROUTINE  BACGRO  ( P * 1 1 . IT  * I AC  * YG * YC * YAD * YU* VC *SDV *NC * SEED  * VU * I REF 
1 *NREF  * VTRAp  * VM  * IM  * 2 ) 

DIMENSION  P(6) *IAC(5) * Y AD ( 5 ) *Y (100) *V(l00) » A ( 100 ) * IR ( 5 ) * X ( 1 00  ) 

1* I TRAP (20) , IREF ( 20 ) * VTRAP(20) 

101  FORMAT  <2X, I3*2X*4(F12.5*2X) ,13) 
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482* 

483* 

484* 

485* 

486* 

487* 

488* 

489* 

490* 

491* 

492* 

493* 

494* 

495* 

49b* 

497* 

498* 

499* 

500* 

501* 

5i'2* 

503* 

5U4* 

5U5* 

80b* 

507* 
508* 
509* 
510* 
511* 
512* 
513* 
514* 
515* 
51b* 
517* 
518* 
519* 
520* 
521* 
5 22* 
523* 
524* 
525* 
526* 
527* 
528* 
529* 
530* 
531* 
532* 
533* 
534* 
535* 
536* 
537* 
538* 


X ( 1 ) =SEED* 1 00000 
NI  = IT-U  + 1 
IR(NCU)=N1 

CALL  RANON  (X#NI>Z>SDV) 

IF  (NC.EQ.o)  GO  TO  5 
DO  4 K=1,NC 

4 IR(K)  = IAC(K)-U  + 1 

5 CONTINUE 

00  6 K=l,NREF  MONDOM 

VTRAP  <K)=0.  MONDOM 

6 ITRAP(K)=IREF(K)-I1+1  MONDOM 

1 = 1 

Y(1)=YC 

IF ( Y ( 1 ) ,LT .P ( 2 ) ) A ( 1 ) =0 . 

IF  (Y(1).GT.P(5>)  A ( 1 ) =P ( 3 )-P(6) * ( Y ( 1 ) “P ( 5 ) ) **2 

IF  ( l Y < 1 ) ,oE.P(2) ) .AND. ( Y ( 1 ) .LE.PI5) ) ) A ( 1 ) =P ( 3 ) *EXP (-P (4 ) * (P ( 5) -Y 
1(1) )**2) 

V(1)=VC 

10=1 

IF  CIPTRD.NE.l)  WRITE  <b»101>  I,Y(1).V(1>»A(1> »YG*I1 

IM=0 

VM  = 0. 


DO  1 1=2  »N1 » 1 
IFINC.EQ.O)  GO  TO  3 
DO  2 K=1,NC 

IF  (II.LE.1R(K)).0R.(I.GT.IR(K+1)))  GO  TO  2 
I0=IR(K) 

YC  = Y ( 1 0 ) 

YG=YAD(K) 

2 CONTINUE 

3 Y(I)=YG+(YC-YG)*EXP( (I0-I)*P(1)) 

IF(Y (I) . LT . P ( 2 ) ) A ( I ) =0 . 

IF  (Y(I ) ,GT.P(5) ) A(I)=P(3)-P(6)*(Y(I)-P(5) )**3 

IF  ( (Y(I ) ,oE.P(2) ) .AND. (Y(I)  .LE.PI5) ) ) A ( I ) =P ( 3) *EXP (-P ( 4 ) * (P (5) -Y 
1(1) )**2) 

V(I)=V(I-1)+.5*(A(I-1)+A(I) ) 

VlJ)=V(I)+x(I) 

IF  (VM.LE.V(I))  IM=I+I1-1 
VM=AMAX1(VM»V(I)) 

I AB= 11*1-1 

IF  (IPTRD.NE.l)  WRITE  (6»101)  I , Y ( I ) . V ( 1 ) , A ( I ) . YG»  I AB 

DO  7 K=1,NHEF  MONDOM 

7  IF  (I.EQ.ITRAP(K)  ) VTRAP(K)=V(I)  MONDOM 

1 CONTINUE 
YU=Y ( N I ) 

VU=V(NI) 

RETURN 


C 

C THE  SUBROUTINE  SPLEND  CARRIES  OUT  THE  CALCULATIONS  ARISING  FROM 

C EITHER  THE  SPLIT  OR  THE  BLEND  AT  WHICH  A FOOD  PACKET  ENDS. 

C 


SUBROUTINE  SPLEND  ( IOP.JU.NU,  YU.VUrWUf  J* Y,  V»W»FR»SEED*HV»HW»NSTOR, 
1MST » LIND ) 

DIMENSION  J(20) »Y (20) »W(20) »V(20) »FR(20)  » X ( 50 ) * HV ( 20 ) * HW ( 20 ) »MST (2 
20) »LIND(20) 
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539* 

101  format  o 

540* 

200  FORMAT  (»  K AN0  LINO (K > * » 10 ( 12 . 2X . 12 *4X ) ) 

541* 

X(1)=SEED*100000 

542* 

N0=2*NU+2 

543* 

CALL  RANDN  ( X *NO*  0 • > 1 . ) 

544* 

60  TO  ( 10  » 20 ) » I0P 

545* 

10  DO  1 K=1,Nu 

546* 

Y(K)=YU 

547* 

IF  (VU.LT.O.)  VU=0. 

548* 

V(K)=VU+HV(K)*A8S(X<K)  ) 

549* 

W(K)=WU*FR(K)*(1.-HW(K)*ABS(X(K+NU)  ) ) 

550* 

1 CONTINUE 

551* 

60  TO  30 

552* 

20  WU=0, 

553* 

XU=0. 

554* 

YU=0. 

555* 

DO  2 K=1,NU 

5bo* 

WU=WU+W  ( K ) 

5b7* 

XU=XU*W(K)*EXP(V(K) ) 

b5b* 

YU=YU+W(K)*Y<K) 

559* 

2 CONTINUE 

560* 

YU=YU/WU 

5&1* 

VU=ALOG(XU/WU) 

562* 

IF  (VU.LT.O.)  VU=0. 

5o3* 

VU=VU+HV(1)*ABS(X(1)  ) 

564* 

WU=WU*(1.-HW<1)*AB5(X(2)  ) ) 

5o5* 

30  CONTINUE 

5o6* 

C 

567* 

c 

FROM  HERE  TO  THE  END  OF  SPeEND  THE  QUANTITY  LlND(K)  IS  FOUND.  THIS 

boo* 

c 

IS  PART  OF  THE  DOMESTIC  MONITOR  SYSTEM.  LlND(K)=l  IF  THE  K-TH 

5o9* 

c 

STORAGE  IS  INVOLVED  IN  THE  SPLIT  OR  BLEND  BEING  DONE  t OTHERWISE 

570* 

c 

LINC(K)=0. 

571* 

c 

572* 

DO  3 Krl.NsTOR 

573* 

3 LINQ(K)=0 

574* 

M1=MST(JU) 

575* 

LINO ( Ml ) r 1 

576* 

DO  4 K-l.Nu 

577* 

J1=J(K) 

576* 

M1=MST(J1) 

579* 

4 LINO ( Ml ) r 1 

560* 

WRITE  (6*200)  (K.LINDtK)  *K=1,NST0R) 

581* 

RETURN 

5o2* 

5b3* 

c 

564* 

c 

THE  SUBROUTINE  SNtpF  CARRIES  OUT  THE  DETERMINATION  OF  WHETHER  THE 

5a5* 

c 

COOKS  PERCEIVE  THAT  SOMETHING  IS  WRONG  WITH  THE  FOOD  THEY  ARE 

5a6* 

c 

WORKING  ON.  THIS  IS  PART  OF  THE  DOMESTIC  MONITOR  SYSTEM. 

567* 

c 

566* 

(>♦  *********************************************************  ****#jM^## 

5o9* 

subroutine  sniff  (V-*wT.Av,SD*xxrKSM) 

590* 

KSM=0 

591* 

C=AV+SD*XX 

592* 

IF  (VT.GT.C)  KSM=1 

593* 

RETURN 

594* 

END 
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Appendix  V:  Executive  Control  Cards 

In  describing  the  control  cards  needed  to  run  this  simulation  on  the  Natick 
Laboratories  UNIVAC  1106  system,  we  shall  begin  by  assuming  the  following: 

(i)  The  program  BACTRAC  is  stored  in  a file  named  MIMOTAPE  on  magnetic 
tape  reel  number  0513. 

(ii)  The  user  wishes  to  enter  a new  set  of  input  data,  which  is  now  stored 
on  cards  and  which  will  become  the  element  ELNAME  after  it  is  stored  in  the  memory. 

Then  the  arrangment  of  control  cards  is  as  follows,  where  we  are  using  0 to  stand 
for  the  number  zero: 

Control  Card  1:  @ RUN  NLMIMO,  GSD030T,  FTLEESIM 

Control  Card  2:  @ ASG,T  MIMOTAPE.,  6C,  0513W 

Control  Card  3:  @ ASG,PU  MIM. 

Control  Card  4:  @ COPIN  MIMOTAPE.,  MIM. 

Control  Card  5:  @ ELT,ILS  MIM. ELNAME 

The  data  deck  follows  control  card  5. 

Control  Card  6:  @ FOR,  US  MIM. BACTRAC,  TPFS.MOD 

Control  Card  7:  @ XQT 

Control  Card  8:  @ ADD,  PE  MIM. ELNAME 

Control  Card  9:  @ PACK  MIM. 

Control  Card  10:  @ PRT,T  MIM. 

Control  Card  11:  @ FIN 
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On  these  and  subsequent  control  cards,  ELNAME  and  MIM  can  be  replaced  wherever 
they  occur  by  any  other  descriptive  names,  having  twelve  or  fewer  letters,  that  the  user 
prefers. 

The  foregoing  control  cards  will  cause  the  input  deck  to  be  listed  on  the  printer, 
then  used  as  the  data  for  the  BACTRAC  program.  This  data  deck  (i.e.  ELNAME)  is 
added  to  the  contents  of  the  file  MIM.  and  retained  in  the  computer  memory. 

The  most  common  situation  is  that  the  user  will  need  to  make  several  computer 
runs,  either  because  he  wants  to  see  the  simulation  results  under  different  conditions  or 
because  the  input  data  deck  has  errors  in  it  that  have  to  be  corrected.  The  arrangement 
of  the  control  cards  on  subsequent  runs  is  similar  to  that  of  the  first  run  (described  above) 
with  the  following  changes: 

(a)  Cards  2 and  4 are  removed 

(b)  Card  3 is  replaced  by 

@ ASG,A  MIM. 

(c)  On  card  5 replace  I by  U in  column  6 

(d)  Immediately  following  card  5,  instead  of  the  data  cards  are  the  correction 
cards,  which  the  user  must  prepare  following  the  instructions  in  the  UN  I VAC  1100 
Operating  System  Manual,  Section  9.5.  These  cards  cause  the  fcomputer  to  make  changes 
in  the  information  stored  in  ELNAME  and  print  its  revised  version.  The  program  is  then 
executed  using  the  revised  data. 

When  the  user  is  finally  finished  using  the  program,  he  must  remove  it  from  the 
computer  memory.  This  is  done  by  means  of  the  following  control  cards 

Control  Card  1:  as  before 

Control  Card  2:  @ ASG,A  MIM. 

Control  Card  3:  @ DELETE  MIM. 
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Control  Card  4: 


@ PRT,T  MIM. 


Control  Card  5:  @ FIN 

If  a user  meets  difficulties  in  running  this  simulation  on  the  Natick  Laboratories 
UNI  VAC  1106  system,  he  should  consult  the  author  of  this  report  or  members  of  the 
Data  Analysis  Office.  In  attempting  to  run  the  simulation  on  a different  computer  system 
the  user  should  first  consult  with  people  who  are  familiar  with  the  system  he  plans  to 
use.  The  executive  control  cards  described  above  will  be  irrelevant  in  that  case. 
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